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Abstract 



The microquasars recently discovered in our Galaxy offer unique opportunity for a deep insight 
into the physical processes in relativistic jets observed in different source populations. We study 
the temporal and spectral evolution of the radio flares detected from the relativistic ejecta in 
the microquasar GRS 1915+105, and propose a model which suggests that these flares are 
due to synchrotron radiation of relativistic electrons suffering adiabatic and escape losses in 
fastly expanding plasmoids (the radio clouds). Analytical solutions to the kinetic equation 
for relativistic electrons in the expanding magnetized clouds are found, and the synchrotron 
radiation of these electrons is calculated. Detailed comparison of the calculated radio fluxes 
with the ones detected from the prominent flare of GRS 1915+105 during March/April 1994 
provides conclusive information on the basic parameters in the ejecta, such as the absolute 
values and temporal evolution of the magnetic field, speed of expansion, the rate of continuous 
injection of relativistic electrons into and their energy-dependent escape from the clouds, etc. 
The data of radio monitoring of the pair of resolved ejecta enable unambiguous determination 
of parameters of the bulk motion of counter ejecta and the degree of asymmetry between them, 
as well as contain important information on the prime energy source for accelerated electrons, 
in particular, may distinguish between the scenarios of bow shock powered and relativistic 
magnetized wind powered plasmoids. Assuming that the electrons in the ejecta might be 
accelerated up to very high energies, we calculate the fluxes of synchrotron radiation up to 
hard X-rays/soft 7-rays, and inverse Compton GeV/TeV 7-rays expected during radio flares, 
and discuss the implications which could follow from either positive detection or flux upper 
limits of these energetic photons. 
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1 Introduction 



The hard X-ray transient GRS 1915+105 discovered in 1992 (Castro-Tirade- et al. 1992) belongs 
to the class of galactic black hole (BH) candidate sources characterized by an estimated mass 
of the compact object M > 3M , strong variability of the electromagnetic radiation, with the 
spectra extending to energies beyond 100 keV (e.g. see Harmon et al. 1994). The invoking 
similarity of these sources with AGNs, the objects powered by the accreting BHs of essentially 
larger mass scales, has recently found a new convincing evidence, after the discovery of apparent 
superluminal jets in GRS 1915+105 (Mirabel & Rodriguez 1994, hereafter MR94), and in 
another hard X-ray transient, GRO J1655-40 (Tingay et al. 1995, Hjellming & Rupen 1995). 
The existence of relativistic jets in these galactic sources, called "microquasars" , makes them 
unique cosmic laboratories for a deeper insight into the complex phenomenon of the jets common 
in those powerful extragalactic objects (Mirabel & Rodriguez 1995). Indeed, being much closer 
to us than AGNs, the microquasars offer an opportunity for monitoring of the jets in a much 
shorter spatial and temporal scales. Furthermore, the parameter L Edd /d 2 , where tL Edd ~ 
1.3 x 1O 38 (M/M ) erg/s is the Eddington luminosity (an indicator of the potential power of 
the source), and d is the distance to the source, is typically by 2-4 orders of magnitude higher for 
the galactic BH candidates than for AGNs. This enables detection of the galactic superluminal 
jets with Doppler factors even as small as S < 1, in contrast to the case of jets in blazars 
(AGNs with jets close to the line of sight), for detection of which strong Doppler boosting of 
the radiation, due to 8 1, is generally needed. Importantly, for the jets with aspect angles 
close to 90° (which is the case in both GRS 1915+105 and GRO J1655-40), the Doppler factors 
of the approaching and receding components are comparable, so both these components may 
be detectable (see Mirabel & Rodriguez 1995). 

In this paper we present the results of theoretical study of the nonthermal flares of GRS 
1915+105, which are observed in the radio (MR94; Mirabel & Rodriguez 1995; Rodriguez et 
al. 1995; Foster et al. 1996) and possibly also infra-red (Sams et al. 1996; Fender et al. 1997) 
bands, and can be expected at higher photon energies as well, if the electrons in the ejecta 
are accelerated up to very high energies (VHE). The main part of the paper is devoted to 
the explanation of the radio flares observed, which provide important information about the 
principal parameters and processes in the radio emitting plasmoids (i.e. the clouds containing 
relativistic electrons and magnetic fields). In Section 2 we deduce, on the basis of qualitative 
consideration, the principal scenario for interpretation of these flares. In Section 3 we obtain 
analytical solutions to the kinetic equation which describes the temporal evolution of relativistic 
electrons in an expanding magnetized medium. In Section 4 we carry out detailed comparison 
of the fluxes of synchrotron radiation emerging from the model radio clouds with the temporal 
and spectral characteristics of the fluxes observed during the prominent March/April 1994 flare 
of GRS 1915+105 (MR94), which provides conclusive constraints for the model parameters of 
the ejecta. The data of radio monitoring of the two (approaching and receding) ejecta present 
particular interest. We show that these data not only enable unambiguous determination of 
the aspect angle and the speed of propagation of the twin ejecta (which perhaps are a little 
asymmetrical, as follows from the observed flux ratio of the counter components), but also 
contain an important information about the prime energy source for continuous energization 
of the plasmoids. 
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For any reasonable assumption about the magnetic field in the ejecta, the energy spectrum 
of relativistic electrons should extend at least up to energies of a few GeV. It is not excluded, 
however, that the electrons in the ejecta are accelerated beyond these energies. Therefore, 
after determination of the range of basic parameters of the expanding plasmoids, in Section 5 
we speculate that in GRS 1915+105, similar to the case of BL Lacs (e.g. Urry & Padovani 
1995; Ghisellini & Maraschi 1997), the electrons may be accelerated beyond TeV energies, 
E > 10 12 eV, and discuss the fluxes of synchrotron radiation which could then extend up to 
X-ray/7-ray energies, as well as the fluxes of the inverse Compton (IC) 7-rays expected in the 
high and very high energies. In Section 6 we summarize the main results of the paper. 

2 Basic estimates 

2.1 Radio fluxes during outbursts 

Observations of GRS 1915+105 show that the radio emission of the source in its active state 
is characterized by the increase of the radio emission at v ~ (1 — 10) GHz from the flux level 
S v < lOmJy to S u ~ 100 mJy (the plateau state, according to Foster et al. 1996), on top of 
which the episodes of radio flares, with durations from days to a month and fluxes up to ~ 1 Jy , 
are superimposed. The spectra in the plateau state are flat or inverted, a r < (approximating 
S v oc z^~ ar ), which is a typical indicative for a strong synchrotron self-absorption in the source 
(e.g. Pacholczyk 1970). The radio flares, which most probably are related with the events of 
ejection of radio clouds, are distinguished by very rapid, less than a day, rise of the fluxes by 
several times, which then decay on timescales of days or more. Remarkably, already at the 
rising stage of the flares the spectra reveal transition from optically thick (a r < 0) to optically 
thin (a r > 0) emission. Initially, when the fluxes are at the maximum level, the spectra are 
typically hard, with a T ~ 0.5, which may then steepen to a r ~ 1. 

Convincing examples of such spectral behavior are given by Rodriguez et al (1995) where the 
detailed data of radio monitoring of GRS 1915+105 at different wavelengths during 5 months 
from 1 December 1993 to 30 April 1994 are summarized. Thus, from Table 1 of Rodriguez et 
al. (1995) it is seen that the spectrum of the flare on December 6.0, 1993, measured with VLA 
radio telescopes simultaneously at 5 different frequencies from 1.47 GHz to 22.49 GHz was 
self-absorbed in the range of 1.47 < v < 4.89 GHz (a T « —0.73), whereas at higher frequencies 
it corresponded to optically thin synchrotron emission, with the mean slope a r ~ 0.52 in the 
range of 4.89 < v < 22.49 GHz. On that day the flux at v « 4.89 GHz was S u « 500 mJy . 
Meanwhile, on November 29.97 the fluxes were at the level of ~ 100 mJy, and the spectrum was 
inverted in the range up to at least 15 GHz. By December 6.6 the fluxes have increased further, 
reaching S„ = 1083 mJy at 3.28 GHz. Comparison with the flux measured simultaneously by 
IRAM telescope at v — 234 GHz (Rodriguez et al. 1995) shows that in this broad range the 
spectrum was 'optically thin', with a T = 0.51, indicating that the source during < 14 hours 
became transparent at 9 cm band. 

An example of gradual steepening of initially hard spectra, with a r ~ 0.5 at the stage of 
flare maximum, is provided by the outburst on 11 December 1993, when the fluxes measured by 
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Nancay telescope at 21 and 9 cm have increased to 1319 and 901 mJy, respectively (cc r ~ 0.45). 
On the next day the fluxes have decreased by factor of 2, while the spectrum steepened to 
a r ~ 0.6 . Gradual steepening of the radio spectra in time is obvious also in the multiwavelength 
VLA data taken on December 14, 16 and 17, when the spectral slope between 4.89 and 22.49 
GHz was changing as a r ~ 0.72, 0.83 and 1.1, respectively (see Rodriguez et al. 1995). 

Similar spectral evolution, although on significantly larger timescales, was observed for the 
prominent strong radio outburst of 19 March 1994, when for the first time a pare of counter 
ejecta were resolved, and superluminal motion in the Galaxy was discovered (MR94). On March 
24 the spectral index at A ~ 3.5 cm was a r ~ 0.5 , which by April 16 steepened to a r ~ 0.8. 
For an estimated distance d = 12.5 ±1.5 kpc to GRS 1915+105 , the apparent velocities of the 
approaching and receding ejecta made 1.25 ± 0.15 and 0.65 ± 0.08 speed of light, respectively, 
using which the real speed (3 ~ 0.92 of the ejecta moving in the opposite directions at an angle 
9 ~ 70° to the line of sight was inferred (MR94). 

Importantly, since both radio emitting components of this strong and long-lived flare were 
resolved, the data of radio monitoring of this flare presented in MR94 contain a unique infor- 
mation about the processes in the jets. Thus, the detection of the receding component enables 
measurement of the flux ratio S a / S T of the approaching and receding ejecta. In the case of 
identical jets, the ratio of the flux densities measured at equal angular separations of the jets 
from the core is 

S r \l-f3cos6j ' 1 j 

where j — 3, if the fluxes are produced in the moving discrete radio clouds, but j = 2 for the 
brightness ratio of continuous stationary jets (e.g. see Lind & Blandford 1985). The radio 
patterns of the March- April flare clearly show motion of two discrete components, but not a 
stationary jet emission. Then, for the deduced values of (3 ~ 0.92 and 9 = 70°, equation (1) 
predicts the flux ratio S a / S T ~ 12, while S a /S T = 8 ± 1 is observed (MR94). This discrepancy 
between the observed and expected flux ratios can be interpreted as an indication of shocks 
propagating with the speed j3 ~ 0.92 in the continuous fluid of radio emitting plasma which 
is steaming with a velocity f3{ ~ 0.7 (Bodo & Ghisellini 1995). However, it is possible also to 
explain (Atoyan & Aharonian 1997) the observed flux ratio in terms of motion of discrete radio 
sources with the speeds coinciding with the deduced pattern speed, if one would allow for the 
twin ejecta be similar, but not necessarily identical, as implied in equation (1). In Section 4 
we will show that the 'discrepancy' between the observed and expected flux ratios are easily 
removed by a small asymmetry between the intrinsic parameters of the ejecta, e.g. by small 
difference in the velocities /3 a and f3 T of the approaching and receding clouds. Then, deviations 
of S a / S T from the 'expected' value can be rather considered as a measure of asymmetry between 
the twin ejecta. 

2.2 Radio electrons 

Generally, gradual transition from self-absorbed to optically thin spectra at progressively lower 
frequencies is a typical indicator of the synchrotron radiation of an expanding radio cloud. 
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Power-law radiation spectra with a r ~ 0.5 imply power-law distribution of relativistic electrons 
AT (7) oc 7"° e with the exponent a c = 2 a r + 1 ~ 2. Since these radio spectra span a broad 
interval of frequencies from ~ 1 GHz up to < 250 GHz, the electron distribution should be 
power-law within some interval of energies 71 < 7 < 72 , with the ratio of Lorentz factors 
72/71 > 10. To estimate the range of basic parameters of radio clouds, such as the magnetic 
field B, size R, speed of expansion f cxp , etc., in this section we will suppose pure power-law 
electron distribution, ^(7) = A j~ 2 , extending from 7x = 1 up to 72 = 7 max , with a sharp cut- 
off above 7 max , i-e. ^(7 > 7 max ) = 0. For this spectrum of electrons A actually corresponds 
to the total number of relativistic electrons, N e = A , and the total energy in the electrons is 
W e = A m c c 2 hi7 max . 

Characteristic frequency of synchrotron photons emitted by an electron with Lorentz factor 
7 in the magnetic field B is v s ~ 0.3 x (1.5z/b7 2 ), where v B = eB sm$/2nm e c, and 1? is the 
pitch angle (e.g. Ginzburg 1979). For the mean sin-# = \j2/3, the frequency 

u s ~ 10 6 5 7 2 Hz , (2) 

where B is in Gauss. Since the magnetic fields expected in the cloud B ~ (0.03 — 0.3) G (see 
below), for production of synchrotron photons at v ~ 300 GHz, electrons with 7 ~ few times 
10 3 are required. In the case of synchrotron origin of the IR jet observed by Sams et al. (1996), 
or of the rapid IR flares detected by Fender et al. (1997) in GRS 1915+105, relativistic electrons 
with 7 ^> 10 4 are to be supposed. However, in the first part of this paper, where we study the 
radio fluxes, we will not speculate with such possibility and will limit the range of relativistic 
electrons under consideration with 7 max = 10 4 . 



2.3 Magnetic fields, size, and speed of expansion of radio clouds 

In the case of an isotropic source with luminosity L'(v') in its rest frame, an observer at a 
distance d would detect the flux S u = 5 3 L'(v/8)/A^d 2 , where 5 = \/l — (3 2 /(l — [3cos6) is the 
Doppler factor of the source moving with a speed (3 at an angle 9 to the line of sight (e.g. Lind 
& Blandford 1985). For distances rf* = <i/12.5kpc ~ 1, the fluxes expected from an optically 
thin cloud are 

S,~2x IO- 39 A 5 1 ' 2 B 3 / 2 v- 1 ' 2 d; 2 mJy . (3) 

Introducing the ratio 77 = w e /% of the energy densities of relativistic electrons w e = 3W c /4:7iR 3 
to the energy density of the magnetic field wb = B 2 /8ir in the cloud, the coefficient A can be 
expressed as A = 2 x 10 50 77 i? 2 i?f 5 , where i?i 5 is the characteristic radius R of the cloud in 
units of 10 15 cm. Since the radio fluxes observed at v ~ 10 GHz reach the level ~ 500 mJy (for 
strong flares), the magnetic field in the emission region should be 

B = 7.7 x 10- 2 5' 1 r 2/7 Ri5 /7 dt /7 S 2 J 7 G , (4) 

where S* = Sioghz/500 mJy. 

The cloud's radius R can be estimated from the requirement of optical transparency of the 
source with respect to synchrotron self absorption. Calculation of the absorption coefficient for 
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the supposed power-law spectrum of electrons is straightforward (e.g. Ginzburg 1979), resulting 
in k u = 4.4 x 10 15 ?7-B 4 ^ 3 cm -1 . Using equation (4), the size R of the cloud can be expressed 
in terms of the opacity r u = R k v at a given frequency v: 



For the ejecta in GRS 1915+105 the Doppler factor 5 ~ 0.6 (MR94), so the optical trans- 
parency t v < 1 at frequencies v ~ (1.4 — 3) GHz requires a magnetic field at the stage of flare 
maximum B ~ (0.2 — 0.3) Gauss, if the energy densities of relativistic electrons and magnetic 
fields in the ejecta would be at the equipartition level rj ~ 1. 

Due to rather weak dependence of B on the opacity r u , the equipartition magnetic field B eq 
could be significantly less than 0.2 Gauss only if r„ < 0.1. However, to reach even r v = 1 in a few 
days after ejection, the clouds should expand with a very high speeds comparable with the speed 
of light. Indeed, from equation (5) follows that the characteristic size of the clouds ejected on 
19.8 March 1994, should increase to R ~ 8 x lO 14 ^ ^ cm by 24.3 March, At ~ 4.5 days after 
the ejection event, when the flare was detected by Nancy telescope (Rodriguez et al. 1995). 
Then, since the intrinsic time in the rest frame of the cloud At' = 5 At, a speed of expansion 
v cxp > 0.1 x c is to be supposed to provide 1+4 GHz < 1. Expansion speeds even higher could 
be needed for the flares which may be weaker, but become transparent during At < 1 day. 
For example, for the radio flare reported by Gerard (1996), apparent evolution of the fluxes, 
measured at 9 and 21 cm, from an inverted spectrum to the one with a T ~ 0.5 occurred between 
9.0 and 10.0 of July 1996, when the flux at 21 cm has increased from Si.4gh z < 30mJy to 
•Si.4GHz = 140 mJy. Using equation (5), an expansion speed v cxp > 0.2c (!) is found. Thus, 
transition of the flares from optically thick to thin spectral forms on timescales of days implies 
that the radio clouds are expanding with velocities f cxp > (0.1 — 0.2) c. This estimate agrees 
with the clouds expansion speed deduced from the observations, and suggests that the lack of 
red-shifted optical lines from GRS 1915+105 might be due to their large Doppler broadening 
(Mirabel et al. 1997). 

2.4 Steepening of the radio spectra: implications 

The observed gradual steepening of the radio spectra from initial a r ~ 0.5 to a T ~ 1 at the stage 
of fading of the flares, requires a steepening of the energy distribution of the parent electrons 
from the power-law exponent a e ~ 2 to a c ~ 3. The steepening of the electron spectral index 
by factor of 1 is usually attributed to the synchrotron (or inverse Compton) energy losses of 
the electrons, -P s (7) oc r ) 2 , and for extragalactic jets the models of this kind were considered in 
a number of papers (e.g., Marscher 1980; Konigl 1981; Bloom & Marscher 1996; etc). In the 
case of GRS 1915+105, however, such an interpretation of the radio data is impossible. Indeed, 
the synchrotron energy loss time of an electron with energy 7 is 




(5) 



(6) 



t s ( 7 ) ~ 7.7 x lO^- 1 ^- 2 s. 



(7) 
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Using equation (2), the time needed for modification of the electron spectra, in the range of 
energies responsible for the synchrotron radiation at frequencies u, is found: 



Thus, in the magnetic fields limited by equation (6), the synchrotron losses of the electrons 
cannot be responsible for the steepening of the radio spectra, since even for frequencies v ~ 
10 GHz they would require timescales larger than few years. 

There are two principle possibilities to explain the observed spectral evolution. The first 
one is to assume that the relativistic electrons are continuously injected into the cloud with a 
spectrum Q(j, t) oc 7~ ain j, with the power-law index a in j = 2 during the first (1 — 3) days, which 
later on steepens to a m j > 3 so as to provide the steepening of the prompt distribution of elec- 
trons iV(7,t) to the power-law index a c ~ 3. The second option is to assume energy- dependent 
escape of electrons from the cloud, while the injection spectrum would not necessarily steepen. 
Note that in this case as well one needs a continuous injection of the electrons, since otherwise 
for any escape-time r esc (7) the escape of electrons would result in a distribution N(j,i) with a 
sharp (exponential) cut-off, instead of power-law modification, above some 7' defined from the 
equation r csc (7') = t (see Section 3). Thus, in both cases a proper interpretation of the spectral 
steepening during the flare implies continuous injection of relativistic electrons into the clouds. 

In principle, neither of those two options can be excluded without specification of <5(t, t) and 
detailed calculation of the electron distribution iV(7,t). However, there are some arguments 
in favor of the escape of relativistic electrons as the main reason for the steepening of radio 
spectra. Indeed, let us suppose that the electron escape is negligible. Then, in order to 
modify the electron distribution from the power-law index a e ~ 2 at t\ ~ (1 — 3) days to 
a e — > 3 at t ~ 2t\ (as observed), the injection of electrons with the steep spectrum a- m j > 3 
should proceed at least with the same rate Q ("/,£) as initially. However, contrary to this 
requirement, strong steepening of the injection spectrum generally indicates that the efficiency 
of particle acceleration is essentially suppressed, so significant decline of the injection rate 
(5(7, t) would have to be expected. Meanwhile, assuming fast escape of relativistic electrons 
from the cloud, the old 'hard spectrum' electrons could be replaced with a freshly injected 
'soft' ones on timescales ~ r csc . Moreover, since iV(7, t) ~ Q(7,£)t csc (for t ^> r csc ), then an 
energy- dependent escape, r csc (7) oc 7 A with A < 1, could provide the required steepening of 
iV(7, t) even for the same hard injection spectrum Q oc 7~ 2 as initially. This is an important 
possibility, since as we shall see later, the decline of S u (i) observed during March- April 1994 
(MR94) implies substantial injection of relativistic electrons into the radio clouds at all stages 
of the flare evolution. 

2.5 Energetics 

The energy Wb of the magnetic field in the radio cloud can be estimated using equation (4): 



The energy of relativistic electrons W e = 77 Hb, therefore the total energy W to t = W e + Wb oc 
^3/7 _|_ ^—4/7^ reac h es the minimum at 77 = 4/3, i.e. around equipartition between magnetic 



t s (u) ~2.4x itfu^B^' 2 s. 



(8) 



Wb ^ 9.9 x 10 41 5- 2 V -' /7 R% 7 d 8 J 7 St /7 erg . 



(9) 
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field and relativistic electrons in the cloud (e.g. Pacholczyk 1970). For the ejecta in GRS 
1915+105 with 5 ~ 0.6, the minimum energy in radio electrons at the stage of flare maximum 
can be estimated as W e ~ 3 x 10 42 erg, which requires continuous injection of the electrons 
with the mean power P inj ~ (1 — 3) x 10 37 erg/s during At' ~ 1 — 3 days. This estimate is in 
agreement with the one obtained by Liang & Li (1995), however, it does not take into account 
the adiabatic energy losses of electrons in the rapidly expanding cloud. As it is shown by 
accurate numerical calculations in Section 4, the adiabatic losses result in an increase of the 
electron injection needed to the level of P in j ~ 10 38 erg/s. Note that the power needed for bulk 
acceleration of the ejecta in GRS 1915+105 is estimated as P je t ~ (1 — 5) x 10 38 erg/s (Meier 
1996). 

The value of P in j ~ 10 38 erg/s corresponds to the minimum power of relativistic electrons to 
be injected into the cloud. In the case of deviation of the magnetic field B in the cloud from its 
'equipartition' value B cq , the energy requirements would be significantly higher. Defining B eq 
as a magnetic field for which 77 = 1, from equation (6) follows that 77 = (P/P cq )~ 17 / 4 . It means 
that a decrease of B by an order of magnitude is to be compensated by an increase of 77 by 
four orders of magnitude, and correspondingly the increase of W e oc rft 7 by almost 2 orders of 
magnitude, to provide the same flux S*. Thus, assuming a magnetic field Bq ~ 0.03 G at times 
t ~ 3 days after ejection (when the cloud's radius R ~ 10 15 cm), instead of B cq ~ (0.2 — 0.3) G, 
one has to suppose P in j ~ 10 40 erg/s, which significantly exceeds the maximum luminosity of 
GRS 1925+105 observed in the X-rays (Greiner et al. 1996; Harmon et al. 1997). So, B > 0.03 
should be considered as a rather conservative lower limit for the characteristic magnetic field 
in the radio clouds at stages of flare maximum. 

Increasing magnetic fields by an order of magnitude beyond equipartition field B cq , we would 
decrease P in j down to 10 36 erg/s. But now the energy in the magnetic field would strongly 
increase, so in order to account for these magnetic fields, we would have to suppose a total 
power in the jet Pj et > 10 40 erg/s. Moreover, for the fields exceeding P eq even by factor of 2, 
the magnetic pressure in the cloud would significantly exceed the one of the relativistic electrons, 
in which case one would have to explain why a strongly magnetized clouds are expanding with 
sub- relativistic speeds (in the opposite case of 77 <C 1 the answer is obvious). 

For the estimated energy W e ~3x lO 42 ^ 3 / 7 erg, and the supposed spectrum ^(7) = A 7~ 2 
in the range 1 < 7 < 10 4 , the total number of relativistic electrons in the cloud is N e = 
Aq ~ 4 x lO 47 ^ 3 / 7 . Assumption of equal amount of protons results in the cloud's mass M c \ ~ 
7 x 10 23 ?7 3 / 7 g. In the equipartition state this estimate is by an order of magnitude smaller than 
the one in Mirabel & Rodriguez (1995). Note, however, that the magnetic field used by Mirabel 
& Rodriguez is B ~ 0.05 G, which is by factor ~ 5 smaller than B cq . For these magnetic fields, 
which cannot be excluded (moreover, seem quite possible, see below), the parameter 77 ~ 10 3 , 
resulting in estimated mass M c \ ~ 10 25 g. 

The kinetic energy of a cloud moving with the bulk Lorentz-factor V ~ 2.5 is estimated 
as Wkin = (r — l)M c ic 2 ~ 10 45 ?7 3 / 7 erg, i.e. essentially more than W e . This is explained 
by the fact that the mean energy per relativistic electron is 7 = ln7 max ~ 10, while the 
assumtion of equal amount of protons in the cloud immediately increases the energy by factor 
of (r — l)m p /m c 7 ~ 300. Then for acceleration of the ejecta in a day or shorter timescales, 
a power Pj Ct ^> 10 40 would be needed. Obviously, this rather uncomfortable requirement 
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significantly softens if one assumes e + — e~ pair plasma (e.g. Liang & Li 1995; Meier 1996). 
Note, however, that it is quite possible to reduce the required power in the jet down to the level 
~ 10 38 erg/s, while still assuming electron-proton plasma in the ejecta. Indeed, if the spectrum 
of the electrons flattens below some 71 < 100 (for B ~ 0.1 G the synchrotron radiation of these 
electrons is below GHz domain), then the total number of electrons and protons may be 2 
orders of magnitude less than estimated above, resulting in the proportional reduction of the 
estimated kinetic energy W^ in of the ejecta. 

2.6 Evolution of the magnetic field in time 

An important information about the evolution of basic parameters in the radio clouds is found 
from the analysis of the observed rate of decline of the flare S v (t) oc during March- April 
1994 (MR94). Generally, S v {t) depends on the temporal evolution of the total number of radio 
electrons and of the magnetic field B(t). Conclusions concerning possible behavior of Bit) 
could be inferred from the condition of 'freezing' of the magnetic field lines into the highly 
conductive fluid. This implies that B/p oc Al, where p is the plasma mass density, and Al is a 
fluid element along the magnetic field line (see Landau & Lifshitz 1963). The plasma density 
depends on time as pit) oc M c \(t) / [R(t)] s . Assuming now that: 

(a) there is no strong turbulent eddies in the cloud, so the length of a fluid element in the 
expanding medium scales as Al(t) oc R(t); and 

(b) the mass of the cloud is constant, 

the magnetic field would be B oc R~ 2 oc t~ 2 , as long as R ~ v exp t. Then from equation 
(3) follows that the radio flux will drop as S u oc N c {t)t~ 3 . Therefore, to account for the 
observed decline of radio flux, one would have to assume an increase of the total number of 
relativistic electrons as fast as N e oc t 1 - 7 . Even neglecting adiabatic energy losses, such a 
behavior of N e implies in-situ acceleration/injection of the electrons with the rate increasing 
as Q(t) oc t ' 7 , resulting in very hard requirements to be imposed on the jet energetics at the 
late stages. Moreover, the dependence of B oc R~ 2 results in Wb oc R 1 , so the total energy of 
the magnetic field would be extremely high, Wb > (10 44 — 10 46 ) erg, if extrapolating backward 
from Wb ~ 10 42 erg at the stage of flare maximum to earlier stages[|, when R -C 10 15 cm. 

To avoid these problems, we have to suppose that at least one of the assumptions (a) or 
(b) above does not hold. If a turbulence is developed in the plasma, then in addition to 
stretching due to cloud's radial expansion, the fluid lines will be stretched in turbulent eddies, 
so Al oc R 1+a with a > 0, resulting in B(R) oc R~ 2+a (a slower decline due to turbulent dynamo 
effect). If the mass of the cloud is not constant, but rather is accumulated in time as M c \ oc t b , 
then B oc t~~ 2+b (additional 'compression' effect). So, the magnetic field can be approximated 
as B oc R~ m oc t~ m , and an index m < 2 would indicate that the magnetic fields are effectively 
created/supplied in/into the radio clouds. Since Wb oc R 3 ~ 2m , a power-law index m < 1.5 will 
eliminate the energetical problems for initial stages of cloud's evolution. On the other hand, in 
the case of m < 1 similar problems appear again, but now for later stages when R ^> 10 15 cm. 

1 We remind that R15 in Eqs. (4) and (5) relates to the size of an optically thin cloud from which unabsorbed 
flux S 1 * is detected, therefore these equations cannot be used for estimation of B at initial stages. 
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Thus, the range of 1 < m < 1.5 seems to be most reasonable for this parameter. Then the 
observed decline S„(t) oc t~ 1,3 can be qualitatively explained, assuming that the total number 
of radio electrons in the cloud increases as N c oc t n with an index n ~ 1.5m — 1.3 > 0. Even 
neglecting adiabatic and escape losses of the electrons, one would need continuous injection 
of radio electrons with Q ("/,£) oc t n_1 , which means that the injection may be stationary or 
gradually decreasing, but substantial at all stages of the flare. 

2.7 The principal scenario 

Summarizing this Section, the following scenario for the radio flares in GRS 1915+105 can be 
proposed. At the first stages of the outburst, the radio emitting plasma forming a pair of radio 
clouds is ejected at relativistic speeds from the vicinity of the compact object, presumably, 
a black hole in the binary. While moving at relativistic speeds f3 ~ 0.9 (at least, for the 19 
March 1994 event) in the opposite directions, the twin clouds (with similar, but not necessarily 
identical parameters) are also expanding with a high speeds, v cxp ~ (0.1 — 0.2)c, to reach a 
size R ~ (0.5 — 1) x 10 15 cm and become optically transparent for synchrotron self-absorption 
in a few days after expulsion. At that times the equipartition magnetic field in the clouds 
B eq ~ (0.2 — 0.3) G is expected, but it may be also several times smaller, which implies a strong 
impact of the electrons on the dynamics of the expanding clouds. Relativistic electrons are 
continuously injected into the emission region with a spectrum Q{j, t) ~ 7~ 2 , and characteristic 
power Pi n j ~ 10 38 erg/s or more, depending on the magnetic field B < B cq . These electrons may 
be due to, e.g., in-situ acceleration at the bow shock front ahead of the cloud, or a relativistic 
wind of magnetized plasmas propagating in the jet region, and pushing the clouds forward 
against the ram pressure of the external medium. An important point is that, simultaneously 
with injection, the electrons also escape from the cloud on energy-dependent timescales t csc (7), 
which modifies the electron distribution N(j, t), and explains the steepening of the radio spectra 
of the fading flare to a r ~ I. The decline of the flare is due to combination of decreasing 
magnetic field in the expanding cloud, decline of the injection rate Q(j,t), as well as adiabatic 
energy losses of radio electrons and their escape from the cloud. 

3 Relativistic electrons in an expanding medium 
3.1 Kinetic equation 

For calculations of the synchrotron and inverse Compton fluxes expected in the framework 
of above described scenario, we have to find energy distribution function N = iV(7, t) of 
the electrons in this essentially nonstationary conditions. The kinetic equation describing the 
evolution of relativistic electrons which we consider in this section, is a well known partial 
differential equation (e.g. Ginzburg & Sirovatskii 1964): 



ON 
~di~ 



$7 



d 



(PN) + Q . 



(10) 
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The Green's function solution to this equation in the case of time-independent energy losses and 
constant escape time t(j, t) = const was found by Syrovatskii (1959). However, in an expanding 
magnetized cloud under consideration we have to suppose that all parameters depend on both 
energy 7 and time t, i.e. Q = Q(j,t) is the injection spectrum, r = r(7, t) is characteristic 
escape time of a particle from the source, and P = P("f : t) = —(d^f/dt) is the energy loss rate. 

Strictly speaking, equation (10) corresponds to a spatially homogeneous source where the 
energy gain due to in-situ acceleration of particles is absent. Actually, however, it has much 
wider applications, and in particular, it seems to be quite sufficient in our case. Indeed, in 
a general form the equation describing evolution of the local (i.e. at the point r) energy 
distribution function / = /(7,r,t) of relativistic particles can be written as (e.g. Ginzburg & 
Syrovatskii 1964): 

^ = div(D r grad/) - div(uj) + |-(P r /) - |-(6J) + ^(dj) , (11) 

where all parameters depend also on the radius-vector r. First two terms in the right side 
of this equation describe diffusive and convective propagation of particles, the last two terms 
correspond to the acceleration of the particles through the first and second order Fermi mech- 
anisms. If there are internal sources of particle injection and sink (such as production and 
annihilation), then the terms similar to the last two ones in equation (10) should be added as 
well. 

Let us consider a source where the region of effective particle acceleration can be separated 
from the main emission region. It seems to be the case for radio clouds in GRS 1915+105, 
where probable site for in-situ acceleration of electrons is a relatively thin region around either 
the bow shock front formed ahead of the cloud, or possibly a wind termination shock formed 
behind of the cloud, in the contact of the relativistic wind in the jet region with the cloud, as 
discussed below in this paper. Meanwhile the main part of the observed flux should be produced 
in a much larger volume V of the post-shock region in the cloud, since the synchrotron cooling 
time of radio electrons is orders of magnitude larger than the dynamical times of the source. 
Since acceleration efficiency (parameters b r and al T ) should significantly drop outside of the 
shock region, after integration of equation (11) over the volume Vo the last two terms can be 
neglected. 

The integration of the left side of equation (11) results exactly in dN/dt. Integration of 
the two propagation terms in the right side of equation (11) gives the net flux of particles, 
due to diffusion and convection, across the surface of the emission region. Obviously, these 
terms are expressed as the difference between the total numbers of particles injected into and 
escaping from the volume Vo per unit time, so the last two terms of equation (10) are found 
(internal sources and sinks, if present, are also implied). At last, integration of the energy loss 
term in equation (11) is reduced to the relevant term of equation (10), where P corresponds 
to the mean energy loss rate per a particle of energy 7, i.e. P = J P r fd 3 r/ J fd 3 r. Under the 
volume Vo of the cloud we understand the region, filled with relativistic electrons and enhanced 
magnetic field, where the bulk of the observed radiation is produced. 

Thus, equation (10) is quite applicable for the study of sources with ongoing in-situ accel- 
eration, as long as the volume Vq, where the bulk of nonthermal radiation is produced, is much 
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larger than the volume AV^ of the region(s) in the source responsible for effective acceleration 
of the electrons. We should mention here that solutions for a large number of particular cases 
of the Fokker-Planck partial differential equation (which includes the term oc d r for stochastic 
acceleration), corresponding to different combinations of terms responsible for time-dependent 
adiabatic and synchrotron energy losses, stochastic and regular acceleration, have been long ago 
obtained and qualitatively discussed by Kardashev (1962). However, only solutions for energy- 
independent escape of relativistic particles were considered, while in our study it is a key feature 
for proper description of the spectral evolution of the radio flares. Another important point 
is that the Fokker-Planck equation generally may contain singularity, so transition from the 
solutions of that equation (if known), which are mostly expressed through special functions, to 
the case of d r — > may not always be straightforward (for comprehensive discussion of the prob- 
lems related with singularities in the Fokker-Planck equation, as well as general solutions for 
time-independent parameters see Park & Petrosian 1995, and references therein). Meanwhile, 
substitution of the acceleration terms by effective injection terms in the regions responsible for 
the bulk of nonthermal radiation, allows us to disentangle the problems of acceleration and 
emission of the electrons, and enables analytical solutions to the first order equation (10) which 
are convenient both for further qualitative analysis and numerical calculations. 



3.2 Time-independent energy losses 

Suppose first that the escape time is given as r = r(7,t), but the energy losses are time- 
independent, P = P("f). The Green's function solution G(7, t,t ) to equation (10) for an arbi- 
trary injection spectrum Nq^) of electrons implies ^-functional injection Q('~f,t) = iV (7) S(t — 
to) at some instant to- At times t > t it actually corresponds to the solution for the 
homogeneous part of equation (10), with initial condition G(j,to + 0, to) — ^0(7), while 
(7(7, to — 0, t ) = 0. Then for the function F = PG this equation is reduced to the form 

dF dF F 

— no) 

dt d( n(C,t) ' 1 ) 

if introducing, instead of the energy 7, a new variable 

c= * )s £m (13) 

where 7* is some fixed energy. Formally, £ has a meaning of a time needed for a particle with 
energy 7 to cool down to energy 7* (so for convenience one may suppose 7* = 1). The function 
ri((,t) = r[e((),t}, where e is the inverse function to (7(7) which expresses the energy through 
C, i.e. 7 = e((). The initial condition for F((,t) reads 

F(C,t )=P[e{Q]N [e{Q] = U(Q. (14) 

Transformation of equation (12) from variables ((, t) to (s — £ + 1, u — i) results in a partial 
differential equation over only one variable for the function Fi(s,u) = F((,t): 

d J± = F ± , (15 ) 

du ti(s — u, u) ' 
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with the initial condition Fi(s,uq) = U(s — uq) found from equation (14). Integration of 
equation (15) is straightforward: 



F 1 (s, u) — U(s — u ) exp 



dui 



do 



n(s - «i,ui) 



(16) 



In order to come back from variables (s, u) to (7, t), it is useful to understand the meaning of 
the function e(s — x) which enters into equation (16) via equation (14) for U and the escape 
function t\ — > r. Since e is the inverse function to g, then for any z in the range of definition 
of this function we have z = g[e(z)}. Then, taking into account that s = ( + t and ( = 5(7), 
for z = s — x we obtain t — x — g[e(s — x)] — 0(7) . For the function g defined by equation (13), 
this equation results in 

rM7,t) d7i 



t 



x 



I 



(17) 



where ^(7, t) corresponds to e(s — x) after its transformation to the variables (7, t). Thus, for 
a particle with energy 7 at an instant t, the function e(s — x) is the energy T x = 1^(7, t) of that 
particle at time x, i.e. it describes the trajectory of individual particles in the energy space. 

Expressing equation (16) in terms of the Green's function G = F/P, the solution to equation 
(10) for an arbitrary r(7, t), but time-independent energy losses of particles is found: 



G(7,Mo) 



^(rto)iVo(r 



to J 



Pin) 



exp 



da; 



to 



t(T x ,x) 



(18) 



Note that this is not a standard Green's function in the sense that the injection spectrum 
was supposed as an arbitrary function of energy ^(7), and not necessarily a delta-function. 
Actually, it describes the evolution of relativistic particles with a given distribution ^0(7) at 
t = to- The solution for an arbitrary continuous injection spectrum is readily found after 
substitution Nq(j) — * Q(j,to)dtQ into equation (18) and integration over dt : 



N(l,t) 



1 



P{l) 



P(T tQ )Q(T t0 ,t )exp 



dx 



to '''(To;; X) 



dt , 



(19) 



with the function T defined via equation (17). In the particular case of time and energy 
independent escape, r(7, t) = const, this solution coincides with the one given in Syrovatskii 
(1959) and Ginzburg & Syrovatskii (1964) in the form of double integral over t and T, if 
equation (19) is integrated over energy with the use of general relations 



dt 



dx 



= P(r x ) , 



(20) 



which follow from equation (17). 

It worths brief discussion of some specific cases of equation (19). Let the escape of particles 
be energy dependent but stationary, r(7, t) — > ^(7), and consider first the evolution of N(j,t) 
when the energy losses are negligible, so T x ~ 7 for any t. Assuming for convenience that the 
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form of injection spectrum does not change in time, i.e. Q("f,t) = Qo{^)q{t) with q(t < 0) = 
(i.e. injection starts at t=0), equation (19) is reduced to 

JV 7, t) = Q 7 )r 7 / # " 7 * e~ 2 dz . (21) 
Jo 

For stationary injection q(t > 0) = 1 the integral results in simple (1 — e~*/ r ), so iV(7,t) ~ 
Qo{l)t until t < t(7), and then the escape of electrons modifies the particle distribution, 
compared with the injection spectrum, as N("y,i) ~ Qo(l) r (l)- 111 the case of t{^) oc 7~ a 
it results in a power-law steepening of the injection spectrum by factor of A. This is a well 
known result of the leaky-box models in the theory of cosmic rays. In the case of nonstationary 
injection, however, the modification of Qi'j) is different. In particular, for an impulsive injection 
q(t) = 5(t) it is reduced to a sharp cut-off of an exponential type above energies 7$ found from 
r(7t) ~ t. We remind, that this qualitative feature has been used in the previous Section to 
argue that for interpretation of the steepening of radio spectra during the flares by escape, one 
needs continuous injection of electrons. 

For a stationary injection of particles, equation (19) can be transformed to the form 



dz 



P z)r z) 



dT (22) 



using equation (20). In the case of r — > oo (absence of escape) and large t, when r = 
To(7,£) — > oo, equation (22) comes to familiar steady state solution for distribution of particles 
in an infinite medium. If the synchrotron (or IC) energy losses of electrons dominate, P = P2I 2 , 
then r = 7/(1 — P2t r ))- In this case N(^, t) ~ tQo(7) until p 2 t / y < 1, and then the radiative 
losses result in a quick steepening of a stationary power-law injection spectrum by a factor 
of 1. Meanwhile, in the case of impulsive injection the modification of the initial spectrum of 
electrons is reduced to an exponential cut-off at 7 > l/p2t (see Kardashev 1962). 



3.3 Expanding cloud 

Energy losses of relativistic electrons in an expanding medium become time-dependent. For 
the radio electrons in GRS 1915+105 only adiabatic losses due to expansion of the cloud are 
important. The adiabatic energy loss rate is given as P a d = v^y/R, where R is the characteristic 
radius of the source, and v is the speed of spherical expansion (e.g. Kardashev 1962). For the 
electrons of higher energies, however, the synchrotron losses may dominate. For the magnetic 
field we suppose B = B (R /R)~ m , where B and R are the magnetic field and the radius of 
the cloud at the instant t . Thus, the total energy losses can be written as 

p = I ( ft + p 4) • < 23 > 

where \i = 2m — 1. For adiabatic losses Pi = v, but we keep the constants p\ and p 2 in the 
parametric form in order to enable other losses with similar dependence on 7 and R as well. 
Here we will suppose that expansion speed v = const, and consider evolution of the particles 
injected impulsively at instant to with the spectrum G(j,to) = Nq(j), as previously. 
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Since the energy losses depend on time via the radius R(t) = Rq + v(t — to), it is convenient 
to pass from variable t to R. Then for the function $ = 7(2(7, R)/R equation (10) reads: 



R 



<9$ 
OR 



7 



d_ 

97 



ai + a 2 



R^J 



- 1 + 



_7_ 

V T 



(24) 



where a\ = pi/v, 0,2 = Vil v i an d for r now we imply the function t{^,R). Transformation of 
this equation from variables (7, R) to (if) = ln(7/i? M ), £ = lni?) results in the equation 
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<9£ dip 



[(^ + ai + 02e^)$i] - 



1 + 



UTi(^,f) 



(25) 



where $1 = $i(^,£) and Ti(ip,£) = r(e^ + ^,e^). The initial condition at £ = hi.Ro reads 
^KV'jCo) — -Ro _1 'No(Roe^) . Thus, we come to the equation formally coinciding with the one 
considered above, with 'time' (£) independent 'energy' (ip) losses P*(ip) — fj, + a± + a,2e^, and 
arbitrary 'escape' function t*(^>, £) = (1 + e^/v Ti) _1 . The solution to this equation is analogous 
to equation (18): 



1 + c*e* 
1 + 



■N (R%e*°)exp 



S e z dz 

a> wri(^ 2 ,z) 



(26) 



where c* = 02/(a* + ^0 = anc ^ ^ = ^x(V')C) is characteristic trajectory of a particle 
in the 'energy' space ip, which is readily calculated from equation (17) for the given P*: 



tfaW', = - M(c* + e-ty^K-*) - c* 



(27) 



Returning now to the variables 7 and R, the evolution of particles with the energy distribution 
iVo(7) at the instant R = R is found: 



G{-y,R,Ro) 



^0 \ 1 ii 



i2 



A^I^Jexp 



R 



dr 



v JRo r(r r ,r) 



(28) 



The energy T r = T r (7,_R) corresponds to the trajectory of a particle with given energy 7 at 
the instant r = i? in the (r, r) plane, and can be represented as T r = 7 A(7, i?, r) where 



A( 7 ,B,r) 



(?)' 



1 + ^ 



(I 



/i+ai 



(29) 



In the formal case of \i + a\ — > 0, when c* = a 2 /(/i + ai) — > 00, equation (29) tends to the limit 
A' = (R/r) ai /[l + a 2 / yR ai \n(R/r)]. Note that for large 7 the radiative losses may limit the 
trajectory of relativistic electrons at some r* > _R , when A(7,_R, r*) — > 00. For these energies 
r Ro in equation (28) should be taken as 00, but £(7, R, R ) = as far as N (oo) = 0. 

In the general case of continuous injection of relativistic particles with the rate Q{^f,t) — > 
Q(7, -R), the evolution of their energy distribution during expansion of the cloud between radii 
R and R > R is found, using equation (28): 



iV( 7 , R) = G( 7 , R,R ) + - ( -= ) A 2 ( 7 , P, r)Q(r n r) exp 



dz 

T(r*,z) 



dr . (30) 
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Here G(7, R, Rq) given by equation (28) describes the evolution of those particles which have 
already been in the cloud by the instant R , while the second term describes the particles being 
injected during the time interval (t ,t). The substitution R = R + v(t — t ) results in an 
explicit expression for iV(7,t). If the source is expanding with a constant velocity v starting 
from t — 0, such a substitution results in formal changes R — > t, Rq — > to, r — > t 1 , and dr = vdt 
in equation (30). If only the adiabatic losses are important, i.e ai — 1 and a2 = 0, equation 
(29) is reduced to a simple A = t/t , and 

NM = i Aro (i 7 ) exp (_£_i£_) + 

ItoZ^^z^' ) P ( / r(t7/x, x)) ^ ^ ^ 

For the energy-independent escape, r(7,t) = r(i), a similar equation can be obtained from 
the relevant Green's function solution found by Kardashev (1962) in his case "stochastic accel- 
eration + adiabatic losses + leakage" (the single combination of terms in that work where the 
escape term was included), if one tends the acceleration parameter to zero. 

It is seen from equation (31) that for a power-injection Q("f) oc 7~ ain i with a m j > 1 the con- 
tribution of the first term quickly decreases, so at t ^> to only contribution due to continuous 
injection is important. This term is easily integrated assuming stationary injection and approx- 
imating r = r 'y~ 5 (t/to) s , with 5 and s > 0. In the case of s < 1, the energy distribution of 
electrons at t 3> r(7, t) comes to N(^y, t) = Q("f) x r(7, t), similar to the case of non-expanding 
source. If s > 1, the condition t ^> r(7, t) can be satisfied only for large enough 7, so only at 
these energies the energy-dependent escape of particles from an expanding cloud can result in 
a steepening of iV(7, t). 

At the end of this section we remind that equation (30) is derived under assumption of 
a constant expansion speed v. However, it can be readily used in the numerical calculations 
for any profile of v(t), if approximating the latter in the form of step functions with different 
mean speeds v\ in the succession of intervals (ti, ti+i) . Similar approach can be implemented, if 
needed, also for modelling of an arbitrary profile of the magnetic field B(R). 

4 Modelling of the March 19 radio flare 

In the framework of the scenario qualitatively described in Section 2, and using the results of 
Section 3, in this Section we study quantitatively the time evolution of the fluxes of 19 March 
1994 radio flare of GRS 1915+105 detected by Mirabel & Rodriguez (1994) and Rodriguez et 
al (1995). This prominent flare presents particular interest from several aspects: 

(a) until now it remains a unique outburst of GRS 1915+105 where both ejecta have been 
clearly resolved in VLA data and parameters of superluminal ejecta calculated (MR94); 

(b) it was very long-lived flare observed during > 40 days; 

(c) it was very strong outburst, with the reported accuracy of flux measurements ~ 5 %. 
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(d) time evolution of the fluxes from both ejecta at ~ 10 GHz, as well as of the total fluxes at 
1.4 — 3.3 GHz are known. 



4.1 Model parameters 

For calculation of the fluxes expected in the observer frame we approximate the principal model 
parameters as the following functions of energy 7 and time t! in the rest frame of spherically 
expanding cloud: 



(i) Speed of expansion is taken as 



v ex Jt') = ^ -r ■ (32) 



For k > this form of v cxp (t'^ enables deceleration of the cloud's expansion at times t' > t cxp . 
The radius of the cloud is 



R{t) ~T^k 




(33) 



For numerical calculations we approximate the profile of v exp (t') with the mean v\ in the suc- 
cession of short time intervals t- 1+ \ — t\ <C t{, as explained above. 

(ii) The mean magnetic field is supposed to be decreasing with the radius R as 

B(t') = B (R/R )- m (34) 

Here R is the cloud's radius at the intrinsic time t' = 2.7 days after ejection (March 19.8). 
This time in the rest frame of the approaching cloud corresponds to the observer's time to = 
t' /5 = 4.8 days when the flare was detected by the VLA telescope (March 24.6). 

(iii) The injection rate of the electrons into each of the twin clouds is taken in the form of 
Q{l,f)= Q (7)?(O>with 

g(t / ) = (l + t'Ainj)- p , (35) 
which enables the decline of injection rate at t' > t in j, and 

Qo(7) °c 7~ Qinj e" 7/7c • (36) 

In this section we do not specify the exponential cut-off energy 7 C , assuming only that 7 C > 10 4 
to account for the radio data. The coefficient of proportionality in equation (36) is chosen so 
as to provide the total flux 655 mJy detected at 8.42 GHz on March 24.6 (MR94). 

(iiii) The escape time r is approximated after following considerations. Let A sc is the mean 
scattering length of the electrons in the cloud. Then the diffusion coefficient D ~ A sc c/3, 
and the characteristic escape time can be estimated as r ~ R 2 /2D ~ 1.5R 2 /c\ sc . If the 
diffusion would be close to the Bohm limit, which corresponds to A sc about of the Larmor 
radius r L = m c c 2 7 jeB of particles, then r oc R 2 ~ m ^. In this case, however, the escape time of 
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7 < 10 3 electrons from the radio cloud with Rq ~ 10 15 cm and Bq ~ 0.1 G would exceed 10 6 yr. 
Thus, in order to provide escape times of order of a day, the diffusion of GeV electrons in the 
clouds should proceed many orders of magnitude faster than in the Bohm limit. Assuming that 
the scattering of the electrons is due to plasma turbulent waves with energy density u> t urb, the 
frequency of collisions, c/A sc , would be oc u> t urb- A widely used approximation u> t urb oc B 2 /8n 
results in r oc i? 2 ^~ m \ However, in principle u> tU rb may drop faster than the energy density of 
the magnetic field. Therefore it seems reasonable to approximate A sc as 



\ SC /R = C x (R/Ro) u ( 7 /7*) z 



(37) 



with the index u being a free parameter. Then the relation u ~ 2m — 1 would correspond to 
the rate of decrease iiWb oc while larger values of u could indicate a faster decay of the 
turbulence. For the normalization energy we take 7, = 2 x 10 3 (i.e. E* = lGeV). Since r 
cannot be less than the time of rectilinear escape of electrons from the cloud, we come to 



r(l,R) 



R 

c 




-A' 



(38) 



For estimation of C\, we remind that at times of flare maximum t' ~ to the electron 
distribution (7, t') in the range of 7 ~ 300 should be not modified by the escape, which requires 
r(300, Rq) > t ~ Rq/vq- Meanwhile at t' ~ few times t the steepening of the radio spectrum 
is significant (see MR94), so r[300, R(3t )] < 3t - These two requirements result in a rough 
estimates C\ ~ (0.3 — 1) if t exp ^> to, and C\ ~ (0.1 — 1) if t cxp < to- Therefore, characteristic 
scattering of the radio electrons should occur on the lengthscales A sc ~ (0.01 — 1) R. Since in 
the case of A ~ 1 the energy dependent term in equation (38) quickly decreases with increasing 
energy, for the electrons with 7 ^> 7* the escape time is energy- independent, r ~ R/c, so rather 
unusual spectra of relativistic electrons in the ejecta can be expected. 

Thus, the principal parameters of our model are: v 0, t cxp and k for the speed of expansion, 
Bq and m for the magnetic field, tj n j and p for the electron injection spectrum, C\, u 
and A for the electron escape time. Some of these parameters, as a in j ~ 2 or A ~ 1, can be 
fixed rather well using the arguments discussed in Section 2. To define the range of freedom 
of other parameters, especially of those which describe the time evolution of the flare, detailed 
quantitative calculations are needed. 



4.1.1 Initial speed of expansion 

To show the range of possible variations of v , in Fig. 1 we present the total fluxes of synchrotron 
radiation which could be expected in the observer frame at t = 4.8 days after ejection of the 
radio clouds, calculated for 3 different expansion speeds t>o, and assuming 2 different magnetic 
fields i?o- It is seen that in the case of B = 0.2 G, when magnetic field and relativistic electron 
energy densities are on the equipartition level, wb — w e , the values of t> ~ (0.1 — 0.2) c 
can explain the observed fluxes. Assumption of B = 0.05 G results in a better agreement of 
calculated and observed fluxes for v = 0.1 c (dashed curves), but the values of t> ~ 0.05 c still 
are not acceptable. 
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In order to show that such high expansion speed of radio clouds is a common feature for 
the flares of GRS 1915+105, in Fig. 2 we present evolution of the fluxes calculated for the flare 
observed on July 8-10 (Gerard 1996), when between July 9 and July 10 apparent transition of 
the spectra in 1.4 — 3.3 GHz band from optically thick to thin forms has occured. All curves 
in Fig. 2 are normalized to the flux 94mJy at 3.31 GHz observed on July 10.0 (the time t — 
in Fig. 2 corresponds to July 8.0). By appropriate choice of the time of ejection it is possible 
to explain both the flux 83mJy detected at 3.31 GHz, and the flux upper limit at 1.41 GHz on 
July 9.0. However, the flux 140 mJy detected on July 10 at 1.41 GHz can be explained only 
assuming high speed of expansion of radio clouds. As it follows from Fig. 2, for v = 0.2 c and 
0.15 c there is a good agreement between the observed and calculated fluxes, whereas in the case 
of t>o — 0.1c the calculated flux Si.4gh z — 105 mJy is significantly below of the ±20mJy flux 
error limits for the Nancay measurements (see Rodriguez et al. 1995). In Fig. 2 the magnetic 
field at the time to — 1 day is B = 0.35 G, which corresponds to the equipartition field for this 
flare. Decreasing P down to ~ 0.02 G, when w c /w B > 10 5 (!), it is possible to fit the measured 
fluxes with t> ~ 0.1 c, but not significantly less than that. Thus, even for relatively weak flares 
one has to suppose very high speeds of expansion of radio clouds, vo — (0.1 — 0.2) c. 

4.1.2 Magnetic field 

Accurate numerical calculations confirm the estimates in Section 2 for the equipartition mag- 
netic field P cq ~ (0.2 — 0.3) G at times of flare maximum, as well as the fact of fast increase 
of the ratio rj = w c /% and required initial injection power P in j = m c c 2 / 7<3o(7)d7 °f rela- 
tivistic electrons in the case of B < B eq . For example, the fluxes corresponding to the case of 
vq = 0.2 c in Fig. 1 require p n j = 1.2 x 10 38 erg/s (in the energy range 7 < 10 4 ) if -Bo = 0.2 G, 
and Pi n j = 10 39 erg/s if B = 0.05 G. Note that for a fixed B , but different values of other 
model parameters the values of p n j may change, within factor of ~ 3, due to different rates of 
adiabatic and escape losses. Thus, for the fluxes shown in Fig. 1 in the case of P = 0.2 G, but 
v = 0.1c, the injection power P in j = 4.9 x 10 37 erg/s is needed. 

In Fig. 3 we show the temporal evolution of the fluxes calculated for a fixed P = 0.2 G 
but 3 different values of the power- law index m in equation (34), and assuming that electron 
escape is negligible. One might conclude from Fig. 3 that for proper interpretation of the 
observed decline of S v (t) the index m is to be about 1.2 — 1.3. However, the upper panel in 
Fig. 3 clearly shows that without energy- dependent escape there is no spectral steepening of the 
fluxes at later stages of the flare, as observed. The dotted curve in Fig. 3, which is calculated 
for m = 1 assuming energy-dependent escape, shows noticeable, but still insufficient, evolution 
of the spectral index a r . Simultaneously, the decline of radio fluxes becomes faster, and now 
well agrees with the observed slope of S v (t). Thus, most probably, the decrease of magnetic 
field in the expanding clouds occurs with the power-law index m ~ 1. Interestingly, for this 
value of m one could expect that during the cloud's expansion the ratio r) = w s /wb remains at 
some fixed level, r](t) ~ const. Indeed, the magnetic field energy density in a cloud, expanding 
with some mean v^, is decreasing oc t~ 2m , while the time behavior of the energy density of 
relativistic electrons can be roughly estimated asw c ocqt/R 3 oc t~ 2 , so rj oc 

f2(m-l) _ const _ 

A slow decline of rj(t) could be expected due to escape of the electrons. Numerical calculations 
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support this conclusion. 



4.1.3 Time profiles of the expansion/injection and escape rates 

Although the dotted curve in Fig. 3 reveals noticeable steepening of the radio spectrum, it does 
not explain the spectral index a r ~ 0.84 measured on April 16 (MR94). Our attempts to fit 
the observed spectra with a constant speed of cloud's expansion, v exp (t) = const, failed. As 
we discussed in Section 3.3, the energy dependent escape of relativistic electrons may result in 
significant steepening of iV(7,t) only on timescales t ^> r(7, t). Meanwhile, the spectral index 
a T ~ 0.5 on March 24 implies t(j, t) > t at t ~ to- Then, for the escape time given by equation 
(38), it is difficult to expect sufficiently fast increase of the ratio t/r on timescales of few times 
t , if the cloud would expand with a constant speed. 

In order to enable a faster escape of the electrons, one has to assume that expansion of 
the clouds at times t > to significantly decelerates, which implies t cxp < t' . Simultaneously, 
the injection rate of relativistic electrons should also decrease, since otherwise due to a slower 
drop of the magnetic field in time (for decelerated expansion) the decline of S u (t) would be 
significantly slower than observed. The study of possible correlations of the model parameters 
in Eqs. (32) and (35) which could explain the evolution of both S u (t) and a T (t), reveals that 
the best fits are reached when t cxp ~ few days, k ~ (0.7 — 1), and more importantly, when 
the time profiles of the injection rate and the speed of expansion are connected with equations 
p ~ 2 k and ti n j ~ (1 — 2) t cxp . In Fig. 4 we show the spectra calculated for t exp = 2 days 
and k = 0.75, assuming different combinations of the ratios t in j/t exp and p/k. Note that for 
both curves corresponding to p = 2k (solid and dotted lines) it is possible to reach a better 
agreement with the observed slope of S u by a better choice of C\ or u. Meanwhile, in the cases 
of p = k and p = 3k (dashed and dot-dashed curves) the discrepancy between the calculated 
and measured fluxes may be reduced only if one assumes t m j tcx P for V — or t; n j 3> tex P 
for p = 3k . Actually, it would correspond to simulation of the relation p ~ 2k by means of 
changing the ratio between t in j and t cxp for q(t) and v cxp (t). 

For the escape rate, the calculations show that generally the value of parameter C\ in 
equation (38) may vary in the range of (0.1 — 1), and the power-law index u should be ~ 1.5, 
which is somewhat larger than 2m — 1 ~ 1. As discussed above, it implies that the turbulent 
energy density in the expanding ejecta decreases faster than wb- 

Thus, the model parameters can be fixed as follows: ~ 2, A « 1, m ~ 1, the expansion 
speed Vq ~ (0.1 — 0.3) c, depending on B ~ (0.03 — 0.3) G at the time t = to, the characteristic 
expansion time t exp ~ few days (in the cloud's rest frame), t in j ~ 1.5 t cxp , and p ~ 2k with 
k ~ (0.7-1). 

4.2 Injection powered by a beam in the jet region ? 

The relation p ~ 2k, with t in j ~ 1.5t cxp , implies an interesting interpretation for the primary 
source of power for continuous injection/acceleration of relativistic electrons. Namely, this 
kind of temporal behavior of q{t') corresponds, in particular, to a scenario when the power 
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Pinj of relativistic electrons injected into the cloud per unit surface area decreases with time as 
Pi n j oc (t')~ 2 . Then at initial stages, when the cloud expands with a constant speed v cxp (t') ~ v , 
the injection rate q(t') oc R 2 Pi n j = const. At times t' ~ t cxp , when the expansion starts to 
decelerate, q(t') also starts to decline. As follows from equation (33), at times t' 3> t C x P the 
radius R{t') oc (t'/t cxp ) 1 ^ k for k < 1, and R(t') oc ln(t' /t cxp ) ~ const for k — 1. Therefore, 
in the case of /c < 1, the total injection power at t' ^> t cxp will be q(t') oc (£' /t cxp )~ 2k , i.e. 
p = 2/c as needed. Interestingly, it is possible to understand also, why i in j in equation (35) is 
somewhat larger than t cxp : an expanding cloud needs some time of order of At' ~ t cxp after 
the beginning of deceleration for a significant deviation of R(f) from the linear behavior, v 1' . 
This interpretation of the time profile of q(t') suggests an injection function in the form 

q h {t>) = R 2 (t')/(v t') 2 . (39) 

Although the discussion above related to k < 1, the time evolution of radio fluxes observed 
during the March/April 1994 flare can be explained with the injection rate given by equation 
(39) even for the values of k somewhat larger than 1. In Fig. 5 we present the fluxes calculated 
for the injection rate %(£') for a fixed t cxp , assuming 3 different values of the index k in the 
range (0.7 — 1.3). The values of v are chosen so as for different k the cloud's radius i? is the 
same. Then in all 3 cases the spectral indexes a r at the instant t practically coincide. 

Remarkably, the decline of the specific injection rate p in j oc t~ 2 suggests a scenario where 
the injection of electrons is powered by a conically beamed continuous flux of energy from the 
central engine. Such a beam may be supposed in the form of relativistic wind of magnetized 
relativistic plasmas and/or electromagnetic waves (Poynting flux) propagating in the conical 
jet region. Then the energy flux density on the surface of a cloud departing from the central 
source at a constant speed (5 would decrease oc t~ 2 , as long as the total energy flux propagating 
in the jet were on a constant level. In this scenario the relativistic electrons may be supposed 
either directly supplied by the beam, or they may be accelerated on the reverse shock front 
terminating the wind, or both of these possibilities. Here we would notice that the most 
probable scenario for explanation of relativistic electrons in the Crab Nebula up to 10 15 eV 
implies their acceleration on the reverse shock, at a distances ~ 0.1 pc, which terminates the 
ultrarelativistic magnetized wind of (predominantly) e + — e~ plasma produced by the pulsar 
(e.g. Kennel & Coroniti 1984; Arons 1996). 

Note that in the framework of the scenario where the bow shock ahead of the ejecta is 
supposed to be responsible in-situ acceleration of radio electrons, interpretation of p lT ^ oc t~ 2 
law is not so straightforward, especially if one takes into account that the ejecta detected during 
the March/April flare of GRS 1915+105 did not show any noticeable deceleration. However, 
we cannot exclude such a possibility, which would require a detailed study, involving a spatial 
profiles of the density and temperature distributions of the ambient gas. Nevertheless, for 
convenience further on we will refer to equation (39) as to the case of "beam injection", and 
will use this injection profile which allows us to reduce the number of free model parameters. 
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4.3 Pair of radio clouds 

So far we considered the temporal evolution of the total fluxes of the March/April 1994 radio 
flare of GRS 1915+105. Below we discuss what kind of information can be found in the 
radio data when both approaching and receding radio components are resolved, as it is the 
case for GRS 1915+10 (MR94), as well as for the second of presently known superluminal jet 
microquasar GRO J1655-40 (Tingay et al. 1995; Hjellming et al. 1995). 



4.3.1 Flux ratio: an asymmetry between ejecta 

When the pair of oppositely moving radio clouds can be resolved, then in addition to the 
equation for the observed angular speed /i a of the approaching cloud, two more equations 
become available, namely, the equation for the angular speed of receding cloud fi r , and the 
equation for the ratio of the flux densities detected from two sources. The first two equations 
are well known, and in general case, allowing for the parameters of the ejecta to be different, 
for a source at a distance d they can be written as 

At sin 9 a c 
1 - (3 a cos 9 a d 

& sin Or c 

^ l + /3 r cos0 r d ' 1 } 

where (3 a>T and 6^ v are the speeds and the angles of propagation of the approaching and receding 
ejecta. Note that for convenience in equation (41) we have substituted 9 r — > 180° — 9 r , so that 
9 a = 6 T for the ejecta moving in strictly opposite directions. 

The equation for the flux ratio of the counter jets being considered, it should be said that 
usually it is supposed in the form of equation (1), where the index j = 3 for the flux ratio 
expected from the pair of discrete radio clouds, while the value of j = 2 gives the brightness 
ratio of continuous stationary jets, both at equal angular separations a = r from the core 
(e.g. Lind & Blandford 1985). Note, however, that equation (1) is valid only if the jets are 
absolutely identical, which implies (3 a = f3 T = (3 and 9 a = 9 r = 9, as well as equal intrinsic 
luminosities L' a = L' r of the jets. In this case Eqs. (1), (40) and (41) are not independent, and 
■predict that (S^/S,-)^ = (/i a / fi T y +olr . Then, for the observed values of /i a = 17 ± 0.4 mas/day, 
/i r = 9.0 + 0.1 mas/day and a r = 0.84 + 0.04 (MR94), in the case of moving radio clouds (j = 3) 
one would expect the flux ratio (S^/SV^ = 13.1 ± 1.7, whereas S a /S r — 8 ± 1 is observed. 
Remaining in the framework of assumption of identical jets, this discrepancy between the 
observed and expected flux ratios of the pair of ejecta in GRS 1915+105 can be explained 
(Bodo & Ghisellini 1995) only if the real speed of radio-emitting plasma (fluid) in the jets 
would be different from the speed (3 ~ 0.92 of the observed radio patterns attributed, in that 
case, to propagation of the shocks in the fluid. 

Meanwhile, it is possible to explain the observed flux ratio in terms of pattern speeds 
coinciding with the speeds of the radio clouds, if we suppose that the twin ejecta are similar, 
but not absolutely identical, allowing for some asymmetry between them (Atoyan & Aharonian 
1997). An asymmetry between the ejecta implies validity of at least one of the inequalities 
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L' a ^ L' r , 9 a ^ 9 T , or j3 a ^ f3 r . Using the relation S(v) = S 3+olr S'(v) between the apparent and 
intrinsic fluxes of a relativistically moving plasmoids, instead of equation (1), we can write the 
equation for the flux ratio of the pair of plasmoids at equal intrinsic times 

SA _S a (T r \ 3+a < fl + p r cos9\ 3+ai L a 



TJ \l-p a cos9l L' r 



r 



where r a r = l/Jl — f3 ai are the Lorentz factors of bulk motion of the ejecta. 



Since a directly measurable quantity is the ratio (S a /S r )<f> at equal angular separations 
a = 0r = 0, using which the flux ratio S a / S T at equal intrinsic times t' a — t' r — t' can be 
deduced, the relation between these two flux ratios is to be found. In the frame of observer 
the ratio of times corresponding to equal angular separation of the sources from the core is 
(ta/trjtj, = fi r /fi a , while the ratio of apparent times t a ^ r = t' '/5 a , r corresponding to the same 
intrinsic time t' is (t a /t r )t' = S r /S a . Using equations (40) and (41), the relation between the 
ratios of the apparent times corresponding to equal angular separations and equal intrinsic 
times is found: 

t a \ . n 



, , Mr) ' (43) 



t' 



where 

r^^sin^, 

r a /5 a sin^ a 1 ; 

Then, for the power-law approximation of temporal evolution of the fluxes oc t~ K , the 
relation between the flux ratios at equal intrinsic times and equal angular separations reads 

^ = a- . (45) 



Or \S T / (f) 

Thus, in the case of symmetrical jets these two ratios coincide, as expected. But generally they 
are somewhat different, depending on the parameters of the asymmetrical ejecta, which can be 
derived from the system of equations (40)-(42) for a given flux ratio at equal intrinsic times 
Sg,/ S T , which then can be compared with the measured ratio (S^/S^V 

Analytical solutions to equations (40)-(42) are found and discussed in detail elsewhere 
(Atoyan & Aharonian 1997). Due to strong dependence of the ratio S a /S T on the ratio of 
Lorentz-factors of the bulk motion r a /r r , equation (42), even the flux ratio S a / S r ~ 6 can be 
easily explained with a very small difference in the speeds of propagation of counter ejecta /3 a 
and f3 T (see Fig. 6). This implies that an asymmetry in the speeds of ejecta can be considered 
as the prime reason for the discrepancy between the measured and 'expected' flux ratios in 
GRS 1915+105, although some asymmetry in the intrinsic radio luminosities L' a and L' T , due 
to somewhat different content of relativistic electrons and/or magnetic fields in the clouds, is 
possible as well. Any significant difference in the angles of propagation 9 a and 9 r is less probable 
(at least, the ejecta move in strictly opposite directions on the sky, see MR94). 

In Fig. 7 we show the time evolution of the fluxes which could be expected from the ap- 
proaching (south) and receding (north) radio clouds, in the case of 3 different flux ratios at equal 
intrinsic times: S a / S r = 7, 9 and 13. It is seen, that although the total fluxes (as previously, 
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normalized to the flux observed on March 24) in all cases coincide, the south-to-north partition 
of these fluxes is essentially different, and for the flux ratio S a / S T = 7 rather good agreement of 
the calculated fluxes with the ones detected from both ejecta is reached. The parameters of the 
bulk motion of approaching and receding ejecta in that case are equal to f3 a = 0.926, f3 T = 0.902 
and 6 = 70.2°, resulting in the Doppler factors 5 a = 0.55 and 5 r = 0.33 . From equation (45) 
follows (for the index k ~ 1.3 deduced from observations, MR94) that the flux ratio S a / S r = 7 
at equal intrinsic times of the clouds corresponds to the flux ratio at equal angular separations 
(S'a/S'r)^ = 8.6, which is in agreement with the measured flux ratio (S^/Sj-)^ — 8 ± 1. 

4.3.2 Synchronous afterimpulses far away from the core? 

It is seen from Fig. 7 that all data points for the total flux of the ejecta, except for the last 
one (on April 30), deviate from the calculated total flux no more than 10%, which is not too 
far from the reported accuracy ~5% of the fluxes measured by the VLA telescopes (MR94). 
The excess ~ 25 % of the total flux on April 30 could be also explained, if one takes into 
account the second ejection event occurred around April 23 (see MR94). Meanwhile, when we 
separate the fluxes measured from the south and north ejecta, the agreement becomes much 
worse. In particular, this relates to the fluxes of both components on April 9 (corresponding to 
t = 21 days), and to the flux of receding component on April 16 (i.e. t = 28 days). Changing 
a little the model parameters, e.g. assuming a faster escape of the electrons, it is possible to 
reach a better agreement with the flux measured from the receding cloud on April 9. However, 
the discrepancy ~ (30 — 40) % for the two other data points cannot be reduced in that way, so 
another explanation of these fluxes is needed. 

We would like to believe that this discrepancy is connected with a significant increase of 
the fluxes between April 4 and 5 detected by the Nancay telescope. It is seen from Table 
1 in Rodriguez et al (1995), that during 24 hours between these 2 days, corresponding to 
t = 16 — 17 days after the ejection event, the fluxes at both 1.41 GHz and 3.3 GHz have suddenly 
increasedby ~ 30%. Although the VLA was not observing GRS 1915+105 on that days, the 
'echo' of that increase could be present in the flux detected by VLA from the approaching 
component on the next observation date, on April 9. The flux from the receding component 
being considered, one should expect a significant delay between the times of observations of that 
event from the south and north radio clouds, if this increase was due to a powerful bidirectional 
afterimpulse from the central source which would have reached the two clouds at equal intrinsic 
times t' a = t' r = t', and therefore, at different apparent times t a = t'/S a and t r = t'/S r . In 
Fig. 8 we show the fluxes calculated for practically the same model parameters as in Fig. 7, but 
assuming that there was an additional short impulse of injection of relativistic electrons into 
both clouds during At' < 1 day around intrinsic time t' = 9 days after ejection (solid curves). 
For the calculated Doppler factors of counter ejecta S a = 0.55 and 5 r = 0.33, this intrinsic 
time corresponds to the apparent times t a = 16.4 days and t r = 27.2 days for the approaching 
and receding clouds, respectively, i.e. just between April 4-5 and April 15-16. It is seen from 
Fig. 8 that the agreement of the calculated fluxes with the measured ones now becomes better. 
Note that the last two data points (on April 30) in Fig. 8 cannot be explained by another 
afterimpulse, since they coincide in time. The excess fluxes on that day are, most probably, 
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connected with the second pair of plasmoids ejected around April 23 (see MR94). 

If such interpretation of the data is not just an artifact, but corresponds to reality, impli- 
cations for the physics of jets may be very important. It would mean that both clouds are 
energized by the central source being even far away from it, at distances 3> 10 16 cm (!), or oth- 
erwise we have to rely upon a mere coincidence of equal intrinsic times, as well as amplitudes, 
of additional injection of relativistic particles from the bow shocks ahead of two counter ejecta, 
which somehow would impulsively increase the rate of transformation of their kinetic energy 
to the accelerated electrons. Energization of the clouds from the central source implies that 
continuous relativistic flux of energy (the beam), in the form of relativistic wind of particles 
and/or electromagnetic fields, should propagate in the jet region. Then the injection of rela- 
tivistic electrons into the clouds could be supplied directly through this wind and/or the wind 
termination reverse shock on the back side of the clouds. 

It should be noted, that in principle the excess of the radio fluxes on April 9 and April 16 
from approaching and receding clouds, respectively, could be explained also by the synchronous 
increase of the magnetic fields in both clouds by ~ 40 % (the dotted curves in Fig. 8), but not of 
the injection rate of relativistic electrons. The principle difference between the 'magnetic' and 
'electronic' afterimpulses consists in different behavior of the spectral index a r : if the injection 
of relativistic electrons proceeds smoothly, the spectral index at a given frequency evolves also 
smoothly, while the 'electronic' afterimpulse can result in significant hardening (variations) of 
the spectral index a T . Thus, the data of multi-frequency radio monitoring of the clouds will be 
able to distinguish between these two options. 

Interestingly, some indication for temporary hardening of the radio spectra during April 
4 to April 6, i.e. coincident with the time of the supposed 'afterimpulse', can be found in 
the data of observations of the March 19 flare at frequencies 1.41 and 3.28 GHz (Rodrtiguez 
et al. 1995). This effect can be seen in Fig. 9 where we compare our calculations with all 
radio data available for the March 19 flare. For the reported accuracy of the Nancy fluxes 
~ 5% (Rodriguez et al. 1990), the agreement between calculated and measured fluxes at low 
frequencies is worse than at 8.4 GHz. At these frequencies some allowance is to be made for 
possible contribution (flickering) of the unresolved central source into the detected total flux, 
which, in particular, may have significant impact on the correct determination of the power- 
law spectral index a r . Nevertheless, it should be said that while at 8.4 GHz the predicted 
fluxes well agree with the VLA data, at low frequencies they agree with the measurements only 
qualitatively, being systematically higher (by ~ 20 %) than the fluxes measured by Nancy. 

Another important feature seen in the Nancy data, is the observed non-smooth decline of 
the flare, with noticeable variation on timescale of days on the S u t plot. If not attributed to 
the flickering of the central source, this may indicate that perhaps the continuous energization 
of the radio clouds is not as smooth as supposed in our model calculations. Obviously, for 
interpretation of the multiwavelength radio data with accuracy ~ 5 %, one needs more accurate 
theoretical models, which would take into account more complex behavior of the injection and 
escape of relativistic electrons, as well as the spatial non-homogeneity of the radio clouds. 
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5 Predictions for nonthermal radiation at high energies 



In the study of radio spectra of GRS 1915+105, in the previous sections we did not specify the 
maximum energies of the accelerated electrons, assuming only that the power-law distribution 
of the electrons injected into the expanding clouds extend to energies beyond several GeV, to 
account for the radio flares of GRS 1915+105 detected. Meanwhile, it is not excluded that 
the relativistic electrons in the jets of microquasars are accelerated to much higher energies, 
similar to the case of jets in the X-ray selected BL Lacs, where the ultrarelativistic electrons are 
shown up via the synchrotron X-rays and inverse Compton TeV 7-rays (e.g. Urry & Padovani 
1995; Ghisellini & Maraschi 1997). Therefore, after determination from the radio data the 
range of model parameters of the ejecta in GRS 1915+105, below we consider the fluxes of the 
synchrotron and IC radiations which could be expected at higher photon frequences, provided 
that the spectra of injected electrons would extend beyond TeV energies. 

5.1 The synchrotron radiation beyond radio domain 

Interpretation of the IR (K-band) jet observed from GRS 1915+105 by Sams et al. (1995) in 
terms of synchrotron radiation suggests that the injection spectrum of the electrons extends 
beyond tens of GeV. Besides, the injection of relativistic electrons should continue with rather 
high rate during the first t ~ 10 days or more, which are needed for the ejecta to reach the 
distances ^> 10 16 cm from the core as observed. The requirement of continuous injection is in 
agreement with the conclusions derived from radio data, however, there are also some problems, 
connected with the interpretation of the absolute values of the fluxes of the IR jet. 

The reported flux of the IR jet corresponds to the level of maximum radio fluxes observed 
during the strong flares, extrapolated to the K-band with the power-law index a r ~ 0.5. This 
typically implies relativistic electrons with power law index a c ~ 2. Meanwhile, although there 
were no observations in the radio band at the time of detection of the IR jet, typically at 
t > 10 days after the ejection the radio flares are already at the decreasing stage, when the 
electron distribution in the GeV range is to be significantly modified, to account for the steep 
radio spectra at this stage. This is demonstrated in Fig. 10 where by heavy lines we show the 
temporal evolution of the spectra of electrons in the clouds, calculated so as to provide the flux 
655 mJy at 3.5 cm on to — 4.8 days after ejection, as previously. Note that the flattening of 
the electron spectra above 10 GeV corresponds to the region where electron escape becomes 
energy independent, being defined essentially by their flight time across the cloud. In Fig. 11 
we show the spectra of synchrotron radiation of these electrons at 3 different times t. It is seen 
that at t — 10 days the synchrotron flux in the IR band is by an order of magnitude less than 
the flux detected by Sams et al (1996). 

To explain the observed high IR flux, we note that the radiation spectra shown by solid 
curves in Fig. 11, correspond to the emission of relativistic electrons inside the clouds with the 
mean magnetic field B(t) and radius R(t). Meanwhile, the electrons escaping from the region 
of high magnetic field, i.e. the cloud, do not immediately disappear, but rather should spend 
some time at distances comparable with R, before departing to large distances from the cloud. 
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In Fig. 10 the thin lines correspond to the energy distribution of the electrons N'^) which have 
already left the cloud, but still are in the region < 3R(t). The calculations are done assuming 
that outside of the cloud the electrons propagate with characteristic radial velocity ~ c/3 
spiraling along the magnetic field lines. The thin lines in Fig. 11 show the synchrotron spectra of 
those electrons calculated for the magnetic field otside the cloud B'(t) = 0.5 B(t). Remarkably, 
at later stages of the flare the synchrotron radiation flux due to the escaped electrons may be 
comparable or even exceeding the radiation of the cloud in the IR/optical region, while in the 
radio domain the radiation produced in the cloud dominates at all timescales. Note that at 
t — 10 days the total nonthermal IR flux predicted in Fig. 11 is only by factor of 3 smaller than 
the one observed. 

The results presented in Fig. 11 indicate that the radiation spectra calculated in the frame- 
work of more realistic, spatially inhomogeneous, model for the entire synchrotron emission 
region could be able to explain the steepening of the radio spectra simultaneously with the 
high IR fluxes much better than the simplified homogeneous cloud model does. Another pos- 
sibility to increase the flux of synchrotron photons in the IR without contradiction with the 
observed steepening of the radio spectra, is connected with the assumption of somewhat harder 
injection spectrum of relativistic electrons, with the power law index a in j ~ (1.8 — 1.9). 

If the electrons in the jets of GRS 1915+105 are accelerated beyond TeV energies, the 
contribution of the synchrotron photons may become important also in the X-ray/7-ray domain. 
Most probable origin of the bulk of observed X-rays in this strong X-ray transient source, with 
the peak luminosities during the flares exceeding 10 39 erg/s (e.g. Greiner et al. 1996), is the 
thermal accretion plasma around the black hole. However, since the observed X-ray spectra are 
rather steep, with a typical power-law index a x ~ 3 in the region of tens of keV (Harmon et al. 
1997), the hard synchrotron radiation of the jets may show up at higher photon energies. It is 
seen in Fig. 11 that during the first several hours after ejection, when the cloud is still opaque 
for synchrotron radiation at GHz frequencies, the synchrotron fluxes at energies > 100 keV 
may dominate over the extrapolation of thermal component. This may result in a significant 
flattening of the overall spectrum at E ~ 100 keV. Although an existence of such feature in the 
spectrum of GRS 1915+105 could be seen only marginally (e.g. see Sazonov et al. 1994), in the 
case of the second microquasar, GRO J1655-40, the X-ray spectra clearly extend up to several 
hundreds of keV (see Harmon et al. 1995). Interestingly, hard tails of the X-ray spectra are 
a characteristic feature of other representatives of the population of galactic BH candidates as 
well, such as Cyg X-l (Ling et al. 1997) or 1E1740. 7-2942 (Churazov et al. 1993; Wallyn et al. 
1996). Extending this spectral feature to the extragalactic jet sources, it is worth noting that 
recently variable hard X-ray spectrum well beyond 100 keV has been detected by BeppoSAX 
during the recent strong flare of the BL Lac source Mrk 501 (Pian et al. 1997). 

5.2 Inverse Compton gamma-rays 

Straightforward evidence for acceleration of relativistic electrons beyond TeV energies in the jets 
of GRS 1915+105 could be provided by detection of the IC 7-rays at very high energies (VHE), 
E-y > 100 GeV. Calculations in the framework of synchrotron-self-Compton model show that 
during the strong flares one may expect detectable fluxes of VHE 7-rays, if the magnetic field in 
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the ejected plasmoids would be significantly below of the equipartition level. In Fig. 12 we show 
the spectra of the synchrotron and IC radiations expected from GRS 1915+105 at t — 0.1 day 
after the ejection event, calculated for the same model parameters as in Fig. 11, except for 
somewhat higher exponential cut-off energy, E c = 20 TeV, assuming 3 different magnetic fields 
Bq at the instant to- It is seen that for the magnetic field Bq = 0.2 G, when the ratio of the 
electron to magnetic energy densities, r] = w c /wb, is at the level close to the equipartition 
(see Fig. 13), the fluxes of the IC 7-rays are rather small, whereas assumption of the fields 
B < 0.1 G results in a significant increase of the IC 7-rays . 

Such a strong dependence of the expected IC 7-ray fluxes on the magnetic field is explained 
by strong dependence of the synchrotron radiation flux on the magnetic field, S u oc B 1+ai (e.g., 
Ginzburg 1979). For magnetic fields smaller by factor of 2, one actually requires an increase of 
the injection rate of the electrons by factor of ~ 3 to provide the same (observed) radio flux. 
Therefore, assumption of different magnetic fields B corresponds to assumption of different 
injection power of relativistic electrons, as shown in Fig. 13. The change of the magnetic field 
Bq has an additional strong impact on the intensity of IC 7-rays, since for a given field of the 
soft photons (normalized to the observed radioflux), the increase of the magnetic field by factor 
of a results in the decrease of the ratio of the photon to magnetic field energy densities by 
factor of a 2 . Therefore, the share of the injection power of VHE electrons which is channeled 
into the IC 7-rays , is essentially reduced. 

In Fig. 14 we show the time evolution of the integral fluxes of the IC 7-rays calculated for 
the same parameters as in Fig. 11, but for 2 different energies of the exponential cut-off in the 
spectrum of injected electrons: E c = 20 TeV (solid curves), and E c = 1 TeV (dashed curves). 
The spectra of the hard X-rays/soft 7-rays due to synchrotron radiation of electrons in the case 
of E c = 20 TeV are also shown (dot-dashed curves). The supposed magnetic field B = 0.05 G 
requires the 'beam injection' power of the electrons P ui j ~ 4 x 10 39 erg/s. Fig. 14 shows that in 
the case of acceleration of the electrons in the ejecta beyond 10 TeV, the synchrotron photons 
during the flare dominate up to lOMeV over the flux of IC 7-rays. The flux of VHE 7-rays 
may be on the level of the Crab flux during several hours after ejection of the radio clouds. In 
a few days the flux drops to the level of 0.1 Crab, which is still detectable by current Imaging 
Cherenkov telescopes in the Northern hemisphere (CAT, HEGRA, Whipple), but afterwards 
the source becomes invisible. 

Note that this conclusion does not principally change, even if we take into account the 
contribution due to IC scattering of the electrons after their escape from the cloud (shown in 
Fig. 14 by thin lines, and calculated under the same assumptions as for Fig. 11), which at 
these times may become comparable with the flux produced inside of the clouds. The fast 
drop of the flux of IC 7-rays is connected with the decrease of the total flux of the background 
photons, and fast expansion of the clouds, both resulting in a drastic decrease of the density 
of target photons (most importantly, from the radio to submillimeter domains). Note that 
contribution due to upscattering of the far IR radiation of the dust, which is possibly present 
in the vicinity of GRS 1915+105 (Mirabel et al. 1997), remains very small even if the energy 
density of the dust radiation would be as high as the one of the blackbody radiation, and the 
radiation temperature be small, T ~ 100 K, which is most favorable for the IC production of 
the VHE 7-rays. 
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The fluence which could be expected at high and very high energies on different timescales 
during the strong outburst is plotted in Fig. 15. For the Imaging Cherenkov telescopes, with 
the atmospheric shower collection area > 10 8 cm 2 , one could expect up to several tens of 7- 
rays with E > 300 GeV during few hours of observations (one night) on the first day of the 
outburst. In the range of high energies, E > 100 MeV, the timescales when the fluxes remain 
on a high level, do not exceed several hours. For a detector as EGRET aboard of Compton 
Gamma Ray Observatory, with effective detection area ~ 10 3 cm 2 , the total number of photons 
expected during that time do not exceed few tens. Meanwhile, for the angular resolution of 
the EGRET ~ (3 — 5)°, the number of background photons above 100 MeV due to diffuse 
background (Hunter et al. 1997) in the direction of GRS 1915+105 (I = 45.37°, b = -0.2°) 
during the same time is by factor of 2-3 more. Continuation of the observations beyond the first 
day cannot noticeably improve the statistics of the high-energy photons due to essential drop of 
the signal at t > 1 day after ejection (see top panel in Fig. 15). Thus, the EGRET could detect 
such flares from GRS 1915+105, in the best case, only marginally. However, the future GLAST 
detector, with 2 orders of magnitude higher sensitivity due to essentially larger effective area 
and better angular resolution of the photons (e.g. see Bloom 1996), will be able to detect the 
7-ray flares correlated with (or slightly preceding) the radio flares from GRS 1915+105 even if 
the electrons in the clouds are accelerated only to 10 GeV energy range. 



6 Discussion 

The data of radio monitoring of the galactic microquasars, which provide a unique opportunity 
to follow the evolution of the fluxes from the pair of relativistic ejecta on conveniently short 
time scales, contain large information on the basic parameters and processes in relativistic jets. 

The temporal evolution of the radio flares in GRS 1915+105, with apparent transition of 
the observed fluxes from synchrotron self-absorbed to optically thin spectral forms at the fast 
rising stage, and later on decline of the fluxes on time scales of days, strongly suggests that the 
flares are produced in fastly expanding radio clouds ejected from the core (MR94; Rodriguez et 
al. 1995; Foster et al. 1996). The estimates of the characteristic parameters of the clouds show 
that in a few days after ejection the clouds should reach the size R < 10 15 cm, which implies the 
expansion speeds ~ (0.1 — 0.2) c in the rest frame of the ejecta. This estimate is in agreement 
with the the expansion speed deduced by Mirabel et al. (1997) from the observations of the 
twin radio plasmoids of 19 March 1994 outburst, and implies extremely large broadening of the 
emission lines produced, if any, in the clouds, which could effectively prevent their detection in 
the IR band. 

Estimates of the equipartition magnetic field in the clouds result in the values B cq ~ 0.2 — 
0.3 G at the stage of maximum of radio flares. The synchrotron cooling time of GeV electrons in 
such magnetic fields is orders of magnitude larger than characteristic timescale of days during 
which the radio spectra steepen from the power-law index a r ~ 0.5 to a r ~ 1. Thus, radiative 
losses cannot be responsible for this effect. We argue that for explanation of the observed 
steepening of radio spectra, a continuous injection of relativistic electrons into the radio clouds 
is needed, the prime reason for the fast modification of the energy distribution of the electrons 
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being the energy dependent escape (due to diffusion and/or drift) of the electrons from the 
clouds. Note that simultaneous injection and escape of the electrons implies that these two 
processes proceed from different sides of the cloud, presumably, in the axceal direction for the 
injection and predominantly transverse direction for the escape. 

In the case of magnetic fields close to the equipartition level, B ~ B cq , the energy Wb+Wc1 ~ 
10 43 erg accumulated in the form of relativistic electrons and magnetic fields in the clouds in a 
few days after ejection is needed. This implies continuous injection of relativistic electrons with 
the initial power P in j ~ 10 38 erg/s. However, the injection is to be essentially more powerful, 
-Pinj > 10 39 erg/s, if the magnetic fields were smaller than the equipartition field by a factor of 2 
or more. In that case the pressure of relativistic electrons would essentially dominate over the 
magnetic energy density in the clouds, and could explain, at least qualitatively, fast expansion 
of the clouds with subrelativistic speeds v cxp > 0.1 c. 

Significant information about the processes in relativistic ejecta in GRS 1915+105 becomes 
available after detailed comparison of the radio data of the prominent 19 March 1994 outburst 
(MR94) with the results of accurate model calculations of temporal evolution of the fluxes, 
produced by relativistic electrons in fastly expanding magnetized clouds. The observed rate 
of decline of the radio fluxes can be explained if the magnetic field declines with increasing 
radius of the cloud as B oc Rr m with the index m ~ 1 (see Fig. 3). This implies that the 
total energy of the magnetic field in the cloud is increasing as W-q oc R, so it should be either 
effectively created in the cloud (presumably, due to turbulent dynamo action) or supplied from 
outside. Interestingly, a similar dependence of the magnetic field in the conical jet, B oc r _1 , 
but where r is the jet cross-section radius (and not the cloud's radius, as in our case) have been 
found by Ghisellini, Maraschi & Treves (1985) for the BL Lac objects. Obviously, these two 
approximations result in the same behavior of the cloud's magnetic field on time, B oc 1/t, as 
long as the cloud expands with a constant speed v exp (t) ~ Vq, but for timescales t > t cxp the 
decline of B(t) becomes different (slower). 

The steepening of the radio spectra from the power-law index a T ~ 0.5 on March 24 to 
the index a r ~ 0.84 on April 16 can be explained, if both the injection of the electrons and 
expansion of the cloud would decelerate in a few days after ejection. For interpretation of the 
temporal evolution of both the spectral index and the fluxes of the March 19 outburst, the 
model parameters in equations (32) and (35), which describe the time profiles of the expansion 
speed of the clouds, v exp (t p rime) , and of the injection rate of relativistic electrons, q(t' ), should 
be connected as t in j ~ (1 — 2) t cxp and p ~ 2 k, with k ~ (0.7 — 1) (see Fig. 4). Remarkably, 
relations of this kind can be provided, if we assume that the electron injection is powered 
by the continuous beam of energy, in the form of magnetized relativistic wind of particles 
and/or electromagnetic waves, propagating in the region of conical jet. Fig. 5 shows that the 
'beam injection' profile qh{t) can readily explain the temporal evolution of the total radio fluxes 
observed at 8.42 GHz (MR94). 

The data of radio monitoring of the pair of counter ejecta in GRS 1915+105 not only allow 
determination of both the speed (3 ~ 0.92 and angle 9 ~ 70° of propagation of the ejecta 
(MR94), but also contain principal information of a new quality, which cannot be found in the 
radio data when only the approaching jet is detected. Indeed, interpretation of the measured 
flux ratio (S a /S T )<j, ~ 8 of the pair of radio condensations in terms of a real motion of a pair 
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of radio clouds implies some asymmetry between the intrinsic parameters of the jets in GRS 
1915+105. Fig. 6 shows that, e.g., a small difference in the speeds of propagation of the ejecta 
would be enough to account even for a flux ratio (S' a /S' r )^ ~ 6, which is by factor of 2 smaller 
than the flux ratio to be expected from identical counter ejecta. Asymmetrical ejection of 
the pair of plasmoids generally would require, just from the momentum conservation law, a 
transfer of a significant recoil momentum to the core of ejection, which is the third object in 
the interacting system "two jets + core". Comparison of this momentum with the integrated 
momentum (absolute values) of the gas orbiting in the inner accretion disk, which is the most 
probable site responsible for production of relativistic ejecta (e.g., Blandford & Payne 1982; 
Begelman et al. 1984), suggests (Atoyan & Aharonian 1997) that asymmetrical ejection of pair 
of plasmoids would be able to induce significant structural changes, or even destruction, of the 
inner disk, resulting in a temporary reduction/termination of the fuel supply into this region 
responsible for the thermal X-rays. 

Thus, in the framework of this scenario, the onset of subcritical/supercritical accretion 
would correspond to the active state of the source, with the X-ray flares which can proceed 
both with and without powerful ejection events and, therefore, observable radio flares (Foster 
et al. 1996; Tavani et al. 1996). Powerful ejection of the radio emitting material may be 
accompanied by significant destruction of the inner accretion zone, which may lead to a strong 
decline of the X-ray fluxes simultaneously with production of relativistic ejecta. The time lag, 
up to few days, between the decline of the X-ray fluxes and appearance of strong radio flares, as 
observed in both superluminal microquasars GRS 1915+105 (Foster et al. 1996; Harmon et al. 
1997) and GRO J1655-40 (Tingay et al. 1995; Harmon et al. 1995; Hjellming & Rupen 1995), 
corresponds to the time needed for expansion of the clouds to become optically transparent 
with respect to synchrotron self-absorption. Depending on the time needed for the recovery 
of the inner accretion disk after ejection, the radio flares may appear just in the dips between 
subsequent X-ray flares, as frequently observed (e.g. Harmon et al. 1997). 

The data on the radio flux evolution of the pair of resolved ejecta contain important infor- 
mation on the mechanisms of continuous energization of the plasmoids. The interpretation of 
the ~ 30 % excess of the fluxes detected at 8.42 GHz from the approaching and receding ejecta 
of GRS 1915+105 on April 9 and April 16, respectively, as the result of synchronous short-term 
increase ('afterimpulse') of the injection rate of relativistic electrons into both clouds, as dis- 
cussed in Section 4, implies continuous energization of the clouds by the central source. Then, 
the energy is to be supplied up to distances > 3 x 10 16 cm by bidirectional beam, presumably 
in the form of relativistic wind of particles and electromagnetic fields, emerging from the core. 
Note that the discussed above relations between parameters of the injection rate and cloud's 
expansion, which can be provided by the 'beam injection' profile given by equation (39), can 
be regarded as another indication of injection of a real beam into the clouds. 

Confirmation of this scenario by future multiwavelength radio observations, if GRS 1915+105 
will give us an opportunity to trace another pair of powerful and long lived jets similar to 19 
March 1994 event, would suggest significant softening of the enormous energetical requirements 
rising in the conventional relativistic jet scenario which implies the kinetic energy of the bulk 
motion of the ejecta be the energy reservoir for in-situ acceleration of the electrons on the bow 
shocks formed ahead of the ejecta. Indeed, the total energy of the radio electrons injected into 
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the clouds during first few days with the power p n j ~ 10 39 erg/s is about W- m j ~ 3 x 10 44 erg, 
and this energy is by a factor of 3 larger for all duration of the 19 March flare. If this amount of 
energy were due to kinetic energy of the bulk motion, W^, then the latter would be estimated 
as at least by another factor of 3 higher than W- m ^ since there was not revealed any deceleration 
of the bulk motion of the plasmoids beyond April 1994 (see MR94). Thus, W^ n > 3 x 10 45 erg 
would be needed. Assuming that this amount of energy is transferred to the bulk motion 
'impulsively', during a very short timescale (since relativistic speeds imply that in a hundred 
of seconds the ejecta will be far away from not only the inner accretion disk, but the entire 
binary), one has to suppose enormous power for acceleration of the ejecta, Pj Ct 3> 10 43 erg/s. 

Meanwhile, in the scenario which implies continuous energy supply by the beam, the beam 
power is to be about P in j, i.e. comparable with the super-Eddington luminosities observed 
during the X-ray flares (e.g. Greiner et al. 1996). Moreover, such a beam of relativistic energy 
can push the cloud forward against the ram pressure of the external medium, therefore at initial 
stages the plasmoids may move even with some acceleration, until they would acquire enough 
mass and kinetic energy to continue the flight ballistically. Depending on the parameters of 
the external medium, the stage of significant acceleration of the ejecta could last, presumably, 
up to several days, which may be not enough to see this effect. Note that for GRO J1655-40, 
which is essentially closer to us, the possibility to resolve initial acceleration, perhaps, could 
be better. Remarkably, in the beam model the acceleration of electrons may occur through 
the relativistic wind and/or the wind termination shock on the reverse (i.e. facing the central 
engine) side of the cloud , although a possible contribution due to acceleration on the bow 
shock ahead of the cloud should not be excluded as well. 

Predictions for the fluxes expected beyond radio domain due to synchrotron and inverse 
Comton radiation mechanisms essentially depend on the maximum energy of accelerated elec- 
trons. Assumption of the electrons accelerated beyond 10 GeV would be enough for explanation 
of the IR jet of GRS 1915+105 observed by Sams et al (1996), and strongly variable IR flares 
detected by Fender et al (1997). Electrons with energies ~ 1 TeV produce synchrotron X-rays in 
the keV band, but even in the case of injection of relativistic electrons with P in j ~ 3 x 10 39 erg/s 
(this corresponds to B ~ 0.05 G), the flux of these photons cannot exceed ~ 10 % of the thermal 
X-rays of the accretion disk. 

However, synchrotron radiation may show up in the range of hard X-rays/soft 7-rays, if the 
spectrum of electrons extends to energies > 10 TeV. In that case during first several hours (up 
to 1 day) after the ejection event, the hard synchrotron X-ray fluxes in the range of E > 100 keV 
could significantly exceed the steep fluxes of thermal X-rays. In the case of P in j > 10 39 erg/s 
(i.e. if the magnetic fields in the clouds are by factor of 2-4 smaller than the equipartition field) 
the fluxes on the level J(> 100 keV) > 10~ 3 ph/cm 2 s can be expected, which could be observed 
by the detectors like BATSE, SIGMA, OSSE, or BeppoSAX, currently operating in this energy 
range. Meanwhile, if the magnetic and electron energy densities in the ejecta are about of 
equipartition, so the injection rate of electrons is p n j ~ 10 38 erg/s, then the expected maximum 
fluxes of the synchrotron X-rays will be smaller by one order of magnitude. Remarkably, even 
these fluxes will be possible to detect with the forthcoming INTEGRAL mission intended for 
operation in the energy range 0.1 — lOMeV. The lack of detection of variable fluxes on the 
level of /(> 100 keV) > 10~ 4 ph/cm 2 s during the first day after ejection event will mean that 
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the maximum (cut-off) energy of relativistic electrons in the ejecta is below TeV domain. 

Another, and most straightforward, evidence for TeV electrons in GRS 1915+105 could 
be provided by detection of VHE 7-rays. Calculations in the framework of synchrotron-self- 
Compton model show that during strong flares one may expect detectable fluxes of TeV 7-rays, 
provided that the the injection power of relativistic electrons P in j > 10 39 erg/s (i.e. magnetic 
field in the radio clouds significantly below equipartition level). The time evolution of the 
IC 7-ray fluxes shown in Figs 14 and 15, predicts that during first several hours of a strong 
outburst the 7-ray fluxes above several hundred GeV are on the level of 3 x 10~ n erg/cm 2 s, 
which corresponds to the level of VHE 7-ray flux of the Crab Nebula. In a few days the flux 
drops to the level 0.1 Crab which can be still detected by current Imaging Cherenkov telescopes. 
Further on, however, the source is not detectable. 

Assumption of a lower magnetic field would result in higher fluxes. However, as far as 
a magnetic fields significantly less than B ~ 0.05 G would imply injection power of electrons 
-Pinj > 10 40 erg/s, it is difficult to expect 7-ray fluxes essentially exceeding the ones shown in Figs 
14 and 15. On the other hand, in the case of in-situ acceleration/injection of the electrons with 
Pinj ~ 10 38 erg/s and magnetic fields on the level of equipartition (B ~ (0.2 — 0.3) G at times 
when R ~ 10 15 cm), the IC fluxes dramatically decrease. Therefore either positive detection or 
upper limits of VHE 7-ray fluxes, being combined with hard X-ray observations above 100 keV, 
could provide robust constraints on the magnetic fields and efficiency of acceleration of electrons 
beyond TeV energies. 

It is important to note, however, that due to narrow field of view of the Cherenkov tele- 
scopes and possibility to do observations only during few night hours per day, the detection 
of the episodes of VHE 7-ray emission will be not an easy task. Fortunately, the predicted 
duration of high fluxes of the VHE 7-rays is up to few days (see Fig. 15). Then, provided that 
our interpretation of the observed anticorrelation/delay between strong X-ray and radio flares 
given above is the case, the instant of a powerful ejection event, and therefore the time (i.e. 
subsequent first night) most favorable for the VHE observations, can be determined from a 
strong temporary drop of the thermal X-ray fluxes during the high state of the source. Impor- 
tantly, occurrence of each strong ejection event can be checked retrospectively by appearance 
of a strong radio flare which is 'delayed' from the time of ejection by 1-3 days. 

Information about the IC 7-rays can be obtained also in the range of high energies, E > 
100 MeV. It should be noted, however, that even during the first few hours of the flare, the 
fluxes expected in this energy range do not essentially exceed 10~ 6 ph/ cm 2 s. This implies that in 
the optimistic case, only marginal detection of such 7-ray flare can be expected for the EGRET 
which has effective detection area ~ 10 3 cm 2 and angular resolution of few/several degrees. We 
can predict, however, that the future GLAST instrument, with 2 orders of magnitude higher 
sensitivity than EGRET, will be able to see the 7-ray flares which would precede, up to few 
days, strong radio flares from GRS 1915+105. 
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Figure 1: The spectral fluxes expected from GRS 1915+105 at to = 4.8 days after ejection of a 
pair of relativistic plasmoids, (3 = 0.92, in the case of different magnetic fields Bo at the instant 
t = to and different expansion speeds v exp (t f ) = vq- The curves shown by solid, dashed, and dot- 
dashed lines correspond to vo = 0.2 c, 0.1c, and 0.05 c, respectively. Other model parameters 
in Eqs. (32)-(38) are: m — 1, a in j = 2, i in j = 20 days (>> t' Q ), and C\ — > 0. All curves are 
normalized to the flux 655 mJy at v — 8.42 GHz. The fluxes measured by the Nancay and the 
VLA radio telescopes (Rodriguez et al. 1995) on 24 March 1994 are shown by diamonds and 
full dots, respectively. The 5 % error bars for the reported accuracies of these measurements 
are about of the size of the full dots. Note that the VLA measurements have been done 0.3 days 
after the Nancay ones. Our calculations, however, show that the fluxes during that short time 
interval change only by several per cent. 
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Figure 2: The temporal evolution of the fluxes at 1.41 GHz and 3.31 GHz for the radio flare 
of GRS 1915+105 observed by Nancy telescope between July 8.0 and July 10.0, 1996 (Gerard 
1996), calculated for 3 different expansion speeds: Vq = 0.2c (solid curves), 0.15c (dashed 
curves), and 0.1c (dot-dashed curves). Correspondingly, 3 different ejection times - July 8.8, 
July 8.73, and July 8.6, are supposed to fit the fluxes detected (full dots) at 3.31 GHz on July 
9.0 and July 10.0. The instant t = corresponds to July 8.0. The open dots show the flux 
upper limits on July 8 and July 9. 
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Figure 3: The temporal evolution of the fluxes at 1.41 GHz, 3.28 GHz and 8.42 GHz, and of the 
spectral index at 8.42 GHz, calculated for 3 different values of the index m in equation (34), 
assuming negligible escape of the electrons from the cloud which expands with a constant speed 
i>o = 0.2 c. The curves plotted by solid, dashed, and dot-dashed lines correspond to m — 1, 
m = 1.5, and m = 0.5, respectively. The dotted curves are calculated for m — 1 in the case 
of significant energy-dependent escape of the electrons, with parameters A = 1 and C\ = 0.5 
in the equation (38). Other model parameters are the for all curves: B = 0.2 G, a in j = 2, 
t in j = 20 days and p — 1. The full squares correspond to the fluxes of the 19 March 1994 flare 
measured by Nancy telescope (from Rodriguez et al. 1995), and the full dots show the fluxes 
of the VLA (MR94). t — corresponds to the supposed time of ejection (March 19.8). 
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Figure 4: The evolution of the flux S v and spectral index a T at 8.42 GHz expected in the case of 
different relations between the expansion and escape model parameters in equations (32) and 
(35): p = 2k and t in j = 2t cxp (solid lines), p = k and £ in j = t cxp (dashed lines), p = 3k and 
t in j = 2t cxp (dot-dashed lines), p = 2k and t in j = t exp (dotted lines). For all cases k = 0.75 
and t CX p = 2 days are supposed. Other model parameters are : a- m j = 2, vq = 0.2 c, -Bo = 0.3 G, 
m = 1, C x = 0.35, u = 1.5 and A = 1. 
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Figure 5: The evolution of the fluxes and spectral indexes at 8.42 GHz, calculated for the 'beam 
injection' profile of relativistic electrons, qb(t'), assuming 3 different profiles for the expansion 
speed v cxp (t'), with the index k = 1 (solid curves), k = 0.7 (dashed curves), and k = 1.3 
(dot-dashed curves). In all cases the time t exp = 2.5 days, and the initial expansion speeds v 
are chosen so as to provide the same radius of the clouds, Rq — 7.1 x 10 14 cm, at the instant 
to = 4.8 days: t>o = 0.15 c, 0.135 c, and 0.165 c, for the solid, dashed, and dot-dashed curves, 
respectively. The escape parameters C\ = 0.25, u = 1.5 and A = 1, and the magnetic field 
So = 0.1 G are supposed. 



40 



5 6 7 8 9 10 11 12 13 1415 




89 t— - 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

5 6 7 8 9 10 11 12 13 1415 

Flux ratio S Q /S r 



Figure 6: The speeds (bottom panel) and the Lorentz factors (top panel) of the bulk motion 
of the pair of ejecta calculated for a distance d = 12.5 kpc, and angular velocities /i a and fi r 
observed from GRS 1915+105 (MR94), assuming different flux ratios S a /S r . 
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Figure 7: The fluxes expected from the approaching (south) and receding (north) radio clouds, 
in the case of 3 different flux ratios at equal intrinsic times, S a / S r = 7 (solid curves), 9 (dashed 
curve), and 13 (dot-dashed curve), calculated under assumption of equal intrinsic luminosities, 
L' a = L' r , but allowing for the speeds of the bulk motion of the 'twin' ejecta to be different, 
which in these 3 cases are equal to (/3 a = 0.926, (3 r = 0.902), (/3 a = 0.923, (3 r = 0.908), and 
(Pa — Pr — 0.918). The clouds expansion speed v = 0.14c and escape parameter C\ = 0.3 
are supposed. All other model parameters are equal to the ones as in the case of k — 0.7 in 
Fig. 5. The data points plotted by stars correspond to the total fluxes, and the full dots and full 
triangles show the fluxes detected by Mirabel & Rodriguez (1994) from the approaching and 
receding components of the March 19 radio flare. The open dot and open triangle, corresponding 
to the fluxes 525 mJy and 130 mJy, respectively, are plotted for demonstration of a possible 
sharing of the 655 mJy total flux between southern and northern components which were not 
resolved on March 24 (see MR94). 
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Figure 8: The temporal evolution of the fluxes (bottom panel) expected from the southern and 
northern radio clouds of the March 19 outburst in the case of synchronous impulsive increase 
of the injection rate of relativistic electrons ('electronic afterimpulse', solid lines) or of the 
magnetic field ('magnetic afterimpulse', dotted lines) in both plasmoids during At < 1 day 
around intrinsic time t' = 9 days after ejection event. The top panel shows the evolution 
of the spectral index at 8.42 GHz expected for the southern component. The dashed curves 
correspond to the case of smooth behavior of the electron injection and magnetic field (i.e., 
without any afterimpulse). The 'beam injection', with the same model parameters as in Fig. 7 
(for S a /S T = 7), is supposed, except for k — 1 and C\ = 0.4. The error bars on the bottom 
panel correspond to the 5 % accuracy of the measured fluxes. 
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Figure 9: The same as in Fig.8, including the Nancy data at 1.41 GHz and 3.28 GHz (Rodriguez 
et al. 1995). The dotted curves in this figure show the evolution of the total fluxes and relevant 
spectral indexes. 
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Figure 10: Energy distribution of relativistic electrons inside (heavy lines) and outside (thin 
lines) of the approaching radio cloud at the apparent times t = 0.1 day (solid lines), 1 day 
(dashed lines), and 10 days (dot-dashed lines). Typical model parameters are supposed: = 
2, vq = 0.2 c, t cxp = 3 days, k — 1, A = 1, u — 1.5, and C\ = 0.3. The magnetic field at the 
instant t = 4.8 days is equal to B = 0.05 G, and the exponential cutoff energy for the injected 
electrons E c = 10 TeV. 
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Figure 11: The spectra of synchrotron radiation, produced inside (heavy lines) and outside (thin 
lines) of the clouds, by the electrons shown in Fig. 10, at t — 0.1 day (solid lines), t = 1 day 
(dashed lines), and t = 10 days (dot-dashed lines). The cross and the diamond correspond to 
the level of maximum fluxes observed during the flares of GRS 1915+105 at v = 8.42 GHz and 
234 GHz, respectively (Rodriguez et al. 1995). The star shows the flux of the IR jet (Sams et 
al. 1996), and the hatched region shows the level of the hard X-ray fluxes typically detected 
from the source in the energy range > 20 keV during the X-ray flares (e.g. Harmon et al. 1997). 
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Figure 12: The fluxes of the synchrotron (heavy lines) and IC (thin lines) radiations which 
could be expected from GRS 1915+105 at t — 0.1 day after ejection of a pair of radio clouds, 
calculated in the case of exponential cut-off energy E c = 20 TeV, assuming 3 different values 
for the magnetic field at the instant t = t : B = 0.05 G (solid lines) B = 0.1 G (dashed line), 
and B = 0.2 G (dot-dashed lines). All other model parameters are the same as in Fig. 11. 
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Figure 13: The ratios of the electron to the magnetic field energy densities in the expanding 
plasmoids corresponding to the calculations in Fig. 12. Different values of the magnetic field B 
at the instant to (when the cloud's radius R = Rq ~ 10 15 cm) actually imply different injection 
power of relativistic electrons P in j. For the supposed E c = 20TeV, the 'equipartition' magnetic 
field corresponds to B ~ 0.3 G. 
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Figure 14: The fluxes of IC 7-rays expected at different times t from GRS 1915+105 during 
a strong radio flare in the case of exponential cut-off energy E c = 20 TeV (solid lines) and 
E c = 1 TeV (dashed lines). The fluxes produced inside and outside of the radio clouds (see 
text) are plotted by heavy and thin lines, respectively. The dot-dashed curves show the fluxes 
of synchrotron radiation for E c = 20 TeV. 
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Figure 15: The fluences of the high and very-high energy IC 7-rays expected during strong radio 
flares in GRS 1915+105 in the case of 3 different exponential cut-off energies for the accelerated 
electrons: E c = 20 TeV (solid curves), E c = 1 TeV (dashed curves), and E c = 30 GeV (dot- 
dashed curve). All other model parameters are the same as in Fig. 11. The supposed magnetic 
field Bq = 0.05 G implies an injection of relativistic electrons with the initial (i.e. at t <C t- m j) 
power P in j ~ (3 ± 0.6) x 10 39 erg/s, depending on E c . 
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